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Abstract. We investigate the rotational velocities of 
early-type stars in the Sco OB2 association. We measure 
vsini for 156 established and probable members of the 
association. The measurements are performed with three 
different techniques, which are in increasing order of ex- 
pected vsim: 1) converting the widths of spectral lines di- 
rectly to wsinz, 2) comparing artificially broadened spectra 
of low vsini stars to the target spectrum, 3) comparing the 
He I A4026 line profile to theoretical models. The sample is 
extended with literature data for 47 established members 
of Sco OB2. Analysis of the vsini distributions shows that 
there are no significant differences between the subgroups 
of Sco OB2. We find that members of the binary popu- 
lation of Sco OB2 on the whole rotate more slowly than 
the single stars. In addition, we find that the B7-B9 single 
star members rotate significantly faster than their B0-B6 
counterparts. We test various hypotheses for the distribu- 
tion of vsini in the association. The results show that we 
cannot clearly exclude any form of random distribution 
of the direction and/or magnitude of the intrinsic rota- 
tional velocity vector. We also investigate the effects of 
rotation on colours in the Walraven photometric system. 
We show that positions of B7-B9 single dwarfs above the 
main sequence are a consequence of rotation. This estab- 
lishes the influence of rotation on the Walraven colours, 
due primarily to surface gravity effects. 

Key words: stars: early-type; formation; rotation - open 
clusters and associations: individual: Sco OB2 
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1. Introduction 

An outstanding problem in theories of star formation 
is the so-called angular momentum problem. Molecular 
clouds and cloud cores, even if they spin about their axis 
only once per Galactic rotation, contain much more an- 
gular momentum than the stars that form from them (see 
e.g., Spitzer 1968). Observations suggest that an inter- 
stellar parent cloud must lose at least two to four or- 
ders of magnitude of angular momentum for a relatively 
wide binary system and a single star, respectively, to be- 
come dynamically possible (see Mouschovias 1991). How 
does angular momentum get redistributed when stars form 
and what is the resulting distribution of rotational veloc- 
ities? These questions can best be addressed by studying 
young stellar groups in which the stars are unevolved and 
physically related. It is in these groups that the observed 
distribution of projected rotational velocities (vsini) may 
still reflect the distribution of rotational velocities at the 
time of star formation. Observations of the present usim 
distribution provide constraints on the angular momen- 
tum history of forming stars, and contain information on 
mechanisms that lead to angular momentum redistribu- 
tion during star formation. 

Suggested mechanisms for the redistribution of angu- 
lar momentum, include magnetic braking and the forma- 
tion of disks. Magnetic braking is important during the 
early, diffuse stages of star formation (Mouschovias 1991). 
Stellar disks are important during the proto-stellar and 
pre main-sequence phases (see e.g., Bodenheimer et al. 
1993). Another way of redistributing angular momentum 
during star formation is through the formation of binary 
systems. This mechanism can be studied by relating the 
characteristics of the binary population in stellar groups to 
the vsini distribution, which may also provide information 
about tidal interactions in close binaries. Ultimately, ob- 
servations of the vsini distribution in young stellar groups 
may provide information on the star formation process it- 
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self. For a recent review on the role of rotation in star 
formation, sec Bodcnhcimer et al.(1993). 

Another important aspect of the study of rotational ve- 
locities of stars is the effect of rotation on observed stellar 
parameters. It is well known that photometry and spectral 
classification of stars are affected by rotation, due to sur- 
face gravity effects. The resulting misinterpretation of the 
observed stellar parameters will bias age determinations 
for stellar groups and the derived mass distributions (e.g., 
Maeder 1971). 

OB associations are young and mostly unobscured sites 
of recent star formation and as such are well suited for 
studies of the properties of young groups of stars (for re- 
views see Blaauw 1964, 1991). In this paper we focus on 
the Sco OB2 association. Because of its proximity (the dis- 
tance to Sco OB2 is ~ 145 pc, see de Geus et al. 1989) it is 
easily accessible to proper motion studies and as a result 
extensive membership determinations have been carried 
out for this association. Sco OB2 consists of three well- 
known subgroups; Upper Scorpius (US), Upper Centau- 
rus Lupus (UCL) and Lower Centaurus Crux (LCC). The 
Upper Scorpius subgroup is located near the remnants of 
its parental molecular cloud, and some of its stars are still 
partly obscured by gas and dust. We also study a subgroup 
located southeast of US, although the physical reality of 
this subgroup has never been established. 

The ages of the three main subgroups were deter- 
mined most recently from Walraven photometry by de 
Geus et al.(1989). They found ages of 4-5 Myr, 14-15 Myr 
and 11-12 Myr for US, UCL and LCC, respectively. For 
US the kinematic age (i.e., the age determined by trac- 
ing the proper motions of the members back in time) was 
determined by Blaauw (1978, 1991) to be 5 Myr. 

Studies of rotational velocities in Sco OB2 have in the 
past concentrated mainly on US. Slettebak (1968) deter- 
mined vsini for 82 stars in US and UCL and found that 
on average the stars in US rotate faster than those in 
UCL (174 kms -1 vs. 119 kms -1 , but the rms spreads 
on these numbers are ~ 100 kms -1 ) and that the mean 
vsini for B7-A0 stars is larger than that of correspond- 
ing field stars (212 kms -1 vs. 138 kms -1 ). The B0-B6 
stars were found to rotate somewhat more slowly than 
field stars. Rajamohan (1976) derived rotational veloci- 
ties for 112 members of Sco OB2 and found that for stars 
with My < 0™0 the distributions of rotational velocities 
are similar for US and UCL and resemble those of field 
stars. Stars with My > 0™0, all of which are located in 
US, were found to rotate much faster than corresponding 
field stars. 

In this paper we present new projected rotational ve- 
locities from high-resolution spectra for stars that are es- 
tablished or probable members of the Sco OB2 association. 
We briefly describe the observations and sample selection 
in Sect. 2. In Sect. 3 we present and discuss the different 
techniques used to derive vsini from the observed spec- 
tra. In Sect. 4 we discuss the observed distribution of pro- 



jected rotational velocities. Various hypotheses about the 
true distribution of rotational velocities are tested for the 
known members of Sco OB2. We also discuss the interpre- 
tation of the data in the context of these hypotheses. In 
Sect. 5 we combine photometric data with our rotational 
velocities and confirm that rotation affects the colours of 
the stars in the Walraven photometric system. We discuss 
the effects of the colour changes on age determinations 
and derivations of mass distributions. Finally, in Sect. 6 
we summarize our conclusions and suggest future work. 

2. Observations 

The spectra from which the projected rotational veloc- 
ities are derived were all obtained with the ECHELlc 
+ Electronographic Camera (ECHELEC) spectrograph, 
mounted at the Coude focus of the ESO 1.52m telescope 
at La Silla. It consists of a 31.6 grooves/mm echelle grating 
and a 632 grooves/mm grism as a cross-disperser. The de- 
tector was a CCD (thinned, back-illuminated, RCA) with 
640 x 1024 pixels of 15 x 15 /im 2 , and a read-out noise of 65 
e _ . In order to increase the signal to noise on read-out the 
CCD pixels were binned 2x2 (henceforth any reference to 
CCD-pixels is to the binned pixels). The linear dispersion 
of ECHELEC is 3.1 A/mm at 4000 A as measured for our 
data. This implies a velocity scale of 7 kms -1 per binned 
pixel and a resolving power of 21 500. The spectra cover 
the wavelength region 3800-4070 A which was recorded in 
11-12 spectral orders. The formal (i.e., taking only pho- 
ton and read-out noise into account) signal to noise ratio, 
measured at the top of the blaze profile in each order, 
typically varies from <~ 70 to ~ 300 between the blue and 
the red end of the spectrum. Further details on the obser- 
vations are given in Verschueren et al.(1996), which also 
describes the data reduction procedure in detail. 

The data were collected during 1991-1993 in the 
framework of an ESO Key Programme aimed at the de- 
termination of high-precision radial velocities of early-type 
stars in young stellar groups (Hensberge et al.1990). The 
stars that were observed are established or probable mem- 
bers of Sco OB2. They form a subset of a larger sample 
of stars that was submitted for observation by the HIP- 
PARCOS satellite. The details of the selection of the HIP- 
PARCOS sample can be found in de Zeeuw, Brown & 
Verschueren (1994). For the present study a total of 156 
stars was observed, of which 136 already have a measured 
vsini. We note here that the observations were not specifi- 
cally designed for rotational velocity studies. The aim was 
to gather radial velocities and also to observe stars re- 
peatedly in order to detect possible spectroscopic bina- 
ries. However, the data are also well suited for deriving 
a large and homogeneous set of rotational velocities. Be- 
cause of telescope and instrumental limitations the mag- 
nitude limit of our sample is <~ 7™ 5. As a result 75% of 
the stars in our sample are of spectral type earlier than 
B7 (see also Verschueren et al.1996). 
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Of the observed sample 74 stars are established mem- 
bers, inferred from proper motions, found by Blaauw 
(1946) and Bertiau (1958). For 13 stars membership was 
inferred from photometry by de Geus et al.(1989). For the 
analysis in Sects. 4 and 5 we also added vsini data from 
the literature (Slettebak 1968, and Uesugi & Fukuda 1981) 
for 47 established members which were not observed in our 
programme. Of these additional stars, 36 are proper mo- 
tion members, the remaining ones are photometric mem- 
bers. Our observed sample is divided among spectral type 
as follows: 2 O stars, 93 B0-B3 stars, 22 B4-B6 stars, 25 
B7-B9 stars and 14 A-F stars. The additional data from 
the literature consists of 6 B0-B3 stars, 8 B4-B6 stars, 32 
B7-B9 stars and 1 A star. In addition to stars in Sco OB2 
8 early-type stars were observed which are not related to 
the association. A number of these stars were used in the 
process of deriving vsini (see Sect. 3.2). The data were 
collected during four observing runs. The runs are labeled 
E2, E3, E5 and E7 by Verschueren et al.(1996); we will 
use these labels here as well. 



3. Techniques Used to Determine vsmi 

We have used three different methods to obtain vsmi from 
the ECHELEC spectra for the stars in our sample: 

1. Deriving vsmi from the full-width half-maximum 
(FWHM) of weak metal lines. This method is appli- 
cable for vsini < 80 kms -1 . 

2. Deriving vsini from the comparison between an arti- 
ficially broadened sharp-lined template spectrum and 
the spectrum of a star for which vsmi is to be de- 
termined. This method is potentially applicable for 
50 kms -1 < vsmi < 200 kms -1 . This range depends 
on spectral type (see Sect. 3.2 and Appendix C). 

3. Deriving vsini from the comparison of theoretical line 
profiles for the Hei A4026 line with the observed line 
profile. The theoretical profile is taken from atmo- 
sphere models for rotating early-type stars by Collins 
et al. (1991, hereafter CTC). This method is in princi- 
ple applicable for all velocities but is most useful for ve- 
locities for which method 2 fails, specifically for higher 
vsini. 

Methods 1 and 2 rely on the classical model for rotat- 
ing stars, as presented e.g., in Chapt. 2 of Tassoul (1978). 
The third method relies on theoretical models in which 
the line profiles are calculated ab initio. Recently Collins 
& Truax (1995) discussed the consequences of using the 
classical model rather than the physically more appropri- 
ate ab initio model. The most important point Collins & 
Truax make is to beware of the assumptions that underlie 
the classical model for rotating stars. These are: 

A. The observational aspect of a uniformly rotating star 
may be approximated by a circular disk subject to a 



limb-darkening law applicable to all parts of the stellar 
disk. 

B. The limb-darkening law and coefficient appropriate for 
the radiation in a line are the same as for the contin- 
uum. 

C. The form of a line does not change over the appar- 
ent disk; it is simply uniformly Dopplcr-shifted by the 
radial motion resulting from rotation. 

Assumption A means that the radial velocity is con- 
stant along lines parallel to the meridian plane of the star. 
This assumption breaks down at high rotational velocities 
where the surface of the star starts to deform and can no 
longer be described as a sphere. The effects of the defor- 
mation of the star arc twofold. The lines of constant radial 
velocity are not parallel to the meridian plane of the star 
and the line formation process is influenced. 

The second assumption, B, is flawed for the following 
reason. The limb-darkening law describes the variation of 
the specific emergent intensity from the centre of the stel- 
lar disk to the limb. The limb-darkening laws that are usu- 
ally quoted in the literature are derived for the continuum. 
However, because the core of a line has a lower specific in- 
tensity at the centre of the stellar disk, the total variation 
from centre to limb of the intensity at the wavelength of 
the line-core will be less. Hence the limb-darkening coef- 
ficient will be smaller for the line-core. This effect is very 
important in strong lines, where the limb-darkening coef- 
ficient shows a variation from the continuum through the 
wings to the core. 

Assumption C does not hold at high rotational veloci- 
ties because the distortion of the star introduces variations 
in the ionization-excitation equilibrium over the surface 
of the star. Furthermore, the local form of the line is in- 
fluenced by gravity darkening and shape distortion. This 
means that there is a difference between the line profiles 
of an intrinsically slow rotator and those of a fast rotator 
seen pole-on. However, according to calculations in Collins 
& Truax (1995) the differences are very small and of no 
concern to the determination of vsini. 

Before discussing the three methods in detail we should 
mention the format of the spectra. After the data reduc- 
tion process the spectra are in the form of normalized 
flux as a function of CCD-pixel. The pixel scale has to be 
converted to a wavelength scale and for the radial veloc- 
ity studies a log-wavelength scale is used (see Verschueren 
et al.1996). It turns out that due to the specific geometry 
of the ECHELEC instrument, the CCD-pixel spacing in 
the spectra is almost identical to a spacing in In A. For 
methods 1 and 2 we used spectra that were not yet con- 
verted to InA-space and for which the different echelle 
orders were not merged. Furthermore, the continuum of 
the spectra was not yet globally corrected for instrumen- 
tal effects (such as the echelle blaze) , which means we had 
to use a local continuum for method 1. For method 3 we 
only used the order containing the He I A4026 line, and to 
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facilitate the comparison of the data to the model profiles 
the order was rebinned into InA-space. In the rebinncd 
spectra the spacing is 2.8 x 10~ 5 in InA or 8.4 kms -1 . 

3.1. Deriving vsmi for Sharp-Lined Stars 

If one assumes that spectral lines in a star are broadened 
by rotation only, the width of the lines can be used to 
derive vsmi. In practice spectral lines arc also broadened 
by thermal and instrumental effects and by other mecha- 
nisms, such as the Stark effect (see e.g., Gray 1992, Chapt. 
17). Therefore one should use weak metal lines that are in 
principle only affected by thermal and instrumental broad- 
ening; two effects one can correct for. However, the weak 
metal lines can also be broadened due to blending, the 
possible spectroscopic binary nature of the star, or curve 
of growth effects. The latter is important for some of the 
stars in our sample (see below) . Possible spectroscopic bi- 
naries show obvious double lines or they reveal themselves 
as stars for which vsmi changes from one observing run 
to the other. Neglecting the binary nature leads to over- 
estimates of vsini. The presence of line blends is more or 
less automatically taken care of in the conversion of line- 
widths to vsmi described below. 

At low values of i;sini one has to take into account 
the consequences of assumptions B and C. Because the 
lines we use are intrinsically narrow and weak (the me- 
dian residual flux in the line cores is ~ 94%, the aver- 
age and rms spread being ~ 91% and ~ 8%) the limb- 
darkening parameter will not be very different from the 
continuum limb-darkening parameter. The small wave- 
length range of the ECHELEC spectra also ensures that 
the limb-darkening parameter is nearly constant over the 
entire spectrum. 

The widths of the metal lines are measured by fit- 
ting a Gaussian plus a straight-line (pseudo) continuum to 
the spectral line. As much continuum as possible is used 
and where necessary double Gaussians are fitted. The fit- 
ting process returns the FWHM of the best fitting Gaus- 
sian, which is a measure of the second moment of the line 
profile. The second moments for different limb-darkening 
laws are derived in Appendix A; because of our CCD- 
pixcl-spacing we use the formulas for InA-spacing. Note 
that Collins & Truax (1995) give the half-width at half- 
maximum instead of the FWHM. The formal errors on the 
derived FWHM and amplitude of the best fitting Gaus- 
sian vary between <~ 1% and <~ 40%. These errors are 
influenced by uncertainties in the continuum level and a 
line-depth (at line-centre) of at least 4-5% is required if 
both the FWHM and amplitude of the best fitting Gaus- 
sian are to have errors of less than 5%. 

The broadening of the lines depends on the assumed 
limb-darkening law. We used a linear limb-darkening 
law to relate the line-widths to vsmi. Dfaz-Cordoves 
et al. (1995) give limb-darkening parameters for the lin- 
ear, quadratic and square-root approximations to limb- 









1.15 


_ | i i i i | i i i i | i i i i 






1.14 








\ 








C 

ml. 13 
> 








< 








\ 
S 








X 

3=1.12 








1.11 


- I > > > > l > > > > l > > > > 


, , , 1 - 






0.3 0.4 0.5 O.f 
u 


i 0.7 


Fig. 1. The value of the FWHM of an 


infinitely sharp line 


which has been broadened by rotation as 


a function of the lin- 



ear limb-darkening parameter u. The full- width-half-maximum 
is normalized by Ao (see Appendix A) . 



1 




\ 






1/ 




• 


1/ 




y" 












y / 






y / 






*' / 






■ ' / 






" / 






* / 

m / 






*/ 




1 


, , , , 1 


: , 1 





20 40 60 


80 




vsmi (km/s) 





Fig. 2. Analytically calculated FWHM of spectral lines (solid 
curve) vs. the FHWM found by fitting a Gaussian (dots). The 
results are shown as a function of vsmi for an instrumental 
profile of 1.4 pixels FWHM, a thermal broadening of 1 pixel 
FWHM and u = 0.4. The dotted line shows the result of fitting 
FWHM= y/a 2 + b 2 (vsini) 2 to the dots. The fit is determined 
primarily by the points at low vsmi. The error at high vsini 
values amounts to 4% at most. 

darkening. Their values are based on Kurucz (1991) atmo- 
sphere models. The value of the limb-darkening parameter 
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for the linear law varies from <~ 0.3 to <~ 0.7 for the range 
of spectral types in our sample. We chose a constant value 
of 0.4 (characteristic of spectral type B3-B5, which is in 
the middle of our range) to derive vsini from the line- 
widths. Figure 1 shows the variation of the FWHM as a 
function of the limb-darkening parameter. It is clear from 
the figure that only small errors are made when using a 
single value for u. The variation between the extreme val- 
ues is <~ 4%. 

The measured line-widths also contain contributions 
from instrumental and thermal broadening. These contri- 
butions have to be subtracted before vsini is derived. For 
our spectra, the instrumental broadening amounts to 1.4 
CCD-pixels FWHM for run E5 and 1.6 pixels for the other 
runs. These values were derived from the measured widths 
of thorium lines in the wavelength calibration spectra. The 
only metal lines that experience any significant thermal 
broadening (as compared to the instrumental FWHM) are 
the C, N and O lines. The amount of thermal broadening 
in these lines varies from ~ 0.7-1.4 pixels FWHM. We di- 
vided our spectral range roughly into three temperature 
ranges. For stars with temperatures around T e g = 8000 K 
an average thermal profile of 0.7 pixels FWHM is used. 
For T off = 15 000 K and T off = 25 000 K the values are 
0.9 and 1.2 pixels FWHM. The differences are small, but 
they are important for the stars with the sharpest lines. 
In practice the spectral type ranges that correspond to 
the temperatures above are: A1-B9, B8-B4, and B3-B0, 
respectively. 



Table 1. Values of parameters a 2 and b 2 . 
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8000 K 


3.05 


3.40 x 10~ 2 


2.45 


3.35 x lO- 2 


15 000 K 


3.37 


3.40 x lO" 2 


2.77 


3.39 x 10~ 2 


25 000 K 


4.00 


3.39 x 10~ 2 


3.40 


3.40 x 10~ 2 



The fitting of Gaussians to estimate the line-widths 
introduces a bias because of the non-Gaussian form of 
the rotational broadening function. We have determined 
this bias by making artificially broadened Gaussian lines 
that had intrinsic widths corresponding to instrumen- 
tal plus thermal broadening. The Gauss-fitting routine 
was applied for a range of vsini values and the results 
were represented with a function of the form FWHM= 
yja 2 + 6 2 (vsini) 2 . The parameters a 2 and b 2 are then used 
in the derivation of vsini. Figure 2 shows an example 
of this investigation. The FWHM values found with the 
Gauss-fitting routine are compared to the analytically cal- 
culated values (Appendix A). The values for a 2 and b 2 are 
listed in Table 1 for different stellar temperatures for ob- 
serving runs E2/E3/E7 and E5 separately. 

Finally, to convert the widths of metal lines to vsini 
the following procedure is used. The widths are measured 



for as many lines as possible (typically 10-50, depend- 
ing on vsini and spectral type) and converted to a his- 
togram of vsini values for each star. This histogram will 
be contaminated by line-blends as well as by interstellar 
lines. By selecting the bins where most of the vsini val- 
ues are concentrated one can weed out the blended lines, 
which cause vsini to be overestimated, and interstellar 
lines, which lead to underestimates. This method has been 
used before by Day & Warner (1975). In practice it works 
best for sharp lined stars with a large number of lines. 
One can then easily distinguish a peak in the histogram. 
Examples are shown in Fig. 3 for two stars. One is a sharp- 
lined A1V star (vsini = 7 kms -1 ) for which <~ 130 lines 
were measured. The other is a B2V star rotating with 
vsini w 60 kms -1 , and for which 25 lines were measured. 
Also shown is the FWHM vs. line-strength for all lines 
measured. The lines that correspond to the peaks in the 
histograms are indicated as solid circles. Note that these 
lines define a definite locus and that for these two exam- 
ples there is no increase of the FWHM with line-depth. 
All other lines are either blends, interstellar, or affected 
by cosmic rays. In the case of HD 144218 the value of 
vsini is higher and the peak in the FWHM distribution is 
less well defined. 

Once the most likely range for vsini has been selected 
from the histograms we proceed by checking whether there 
is a significant correlation between line-depth and FWHM. 
Such a correlation is expected if lines are on the flat part of 
the curve of growth, where the FWHM increases with line- 
depth. If a significant correlation is found we correct the 
FWHM of each line by fitting a linear model to the data. 
After this correction we check whether the distribution 
of the FWHM of the lines is significantly skewed towards 
high values. If this is the case there may still be line blends 
in the data. We correct for these blends by repeatedly 
discarding the line with the largest FWHM until the rms 
spread in the data changes by no more than 10%. This 
procedure is applied subject to the condition that no more 
than 30% of the lines are discarded. The two quantitative 
criteria were found by trial and error. After the selection 
procedure is completed vsini is calculated as a weighted 
average of the vsini values given by each line. In order 
to ensure that lines with erroneously low formal errors on 
their FWHM do not enter into this average with too much 
weight, the median error in the FWHM is calculated from 
the data and half of this value is added in quadrature to 
all errors. 

3.2. Deriving vsini by Template Broadening 

For stars with vsini larger than about 50 kms -1 the 
method for deriving vsini described above becomes unreli- 
able. This is because the estimated stellar continuum level 
is already significantly too low due to rotational broaden- 
ing, which implies that the measured FWHM of the lines 
is underestimated. Furthermore, at higher values of vsini 
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the number of lines for which one can measure the FHWM 
is reduced, and it is more difficult to identify and eliminate 
line blends. Therefore, a different technique for obtaining 
vsini is required. 

The technique we consider here is that of broadening of 
spectra of sharp-lined stars and comparing the broadened 
spectrum to the spectrum of the target star for which vsini 
is to be determined. Again, the model that is used when 
broadening a template spectrum is the classical model 
for rotating stars. In this case all three assumptions of 
the classical model play a role. The assumptions on limb- 
darkening are now very relevant because of the inclusion of 
broad and strong hydrogen and helium lines in the deter- 
mination of vsini. The limb-darkening we used was again 
linear limb-darkening with a parameter u = 0.4. The ef- 
fects of varying the limb-darkening are discussed at the 
end of Appendix C. At the highest vsini values the defor- 
mation of the target star, which leads to changes in line- 
shapes over its surface, shows up as a mismatch between 
its spectrum and the broadened template spectrum (see 
also Appendix C). The template stars were taken from our 
own sample and they were required to have a low value 
of vsini (~ 10-20 kms -1 ), which was determined using 
method 1. 

Bearing all the mentioned caveats in mind we proceed 
as follows. We broaden a template spectrum for a num- 
ber of values of vsini over a chosen interval and find the 
value for which the difference between broadened template 
and target is minimal. The search interval is chosen based 
on the literature value for vsini or by visual inspection 
of the spectrum of the target. The comparison between 



Fig. 3. Histograms of i;sini values de- 
rived from measured line-widths and 
the line-widths vs. line-strengths. The 
data for HD 144218 (B2V) are shown 
in (a) and (b). The data for HD 72660 
(A1V) are shown in (c) and (d). In- 
dicated with solid circles are the lines 
corresponding to the peaks in the his- 
tograms. For HD 144218 this is the 
velocity range 45kms _1 • vsini < 
65 kms -1 and for HD 72660 the range 
is 0-10 kms^ 1 (the sharpest lines in HD 
72660 formally give negative vsini val- 
ues, after subtraction of thermal and in- 
strumental broadening, and are either 
interstellar, or affected by cosmic rays, 
or under-sampled). 

the broadened template spectrum and the target spec- 
trum is done for 10 orders in the ECHELEC spectra (i.e., 
orders 139-148, see Verschueren et al.1996). To find the 
best vsini value a figure of merit analogous to a \ 2 is con- 
structed for each order separately. This figure of merit is 
the mean absolute difference, over all pixels, between the 
template and the target divided by the noise in the target 
spectrum. This is essentially an automation of the com- 
parison that one would do by eye. Based on the values of 
the figures of merit a weighted average of the vsini val- 
ues found for each order is calculated. Values of vsini at 
the edge of the search interval are discarded. These values 
usually point to noisy orders, orders with low information 
content or orders where there is a significant mismatch 
between the spectra of the template and the target. In all 
cases the values of vsini found for all orders are inspected 
to check whether the results are consistent. The wrong 
choice for the search interval can cause a large spread in 
the data that is not real. 

The advantage of doing the comparison per spectral or- 
der is that it provides information on vsini independently 
for different parts of the spectrum. This provides an effec- 
tive way of detecting mismatches, either intrinsic or due 
to the data reduction process, between the template and 
target spectra by demanding consistency of vsini values. 
A global comparison only gives a single value of vsini and 
spectral type mismatches are then only revealed by visual 
inspection. Furthermore, from the spread in the data and 
a visual inspection of the results one gets a good indica- 
tion of the precision as well as the accuracy with which 
vsini was determined. Finally, there are often weaker lines 
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Fig. 4. Comparison of the target spec- 
trum of HD 121790 (B2IV-V) with the 
best-fit broadened template spectrum 
of HD 96706 (B2V). The thick line 
is the broadened template spectrum. 
The central wavelength of each order 
is indicated, as well as several spec- 
tral lines. The best value of vsini for 
the target star is 118 kms -1 . Note the 
mismatches between the spectra of the 
stars at several places in the continuum 
as well as in the wings of the hydro- 
gen lines. These mismatches can be due 
to intrinsic spectral type mismatches or 
due to the data reduction process. 



present in the broad wings of the Balmer lines. The influ- 
ence of these lines is not taken into account in the mod- 
els by CTC. These lines can make a difference if they 
are strong enough and template broadening automatically 
takes their contribution into account. 

Spectral type mismatches often result from slight dif- 
ferences in temperature or luminosity class. The tempera- 
ture differences are readily apparent in the strength of the 
helium lines for B-stars and the luminosity differences are 
manifest in the hydrogen lines. Another problem to take 
into account when looking for a template spectrum is a 
possible error in spectral classification of either template 
or target. The only way to get around this problem is to 
try different template spectra and find the best match. 

The information on vs'mi in the rotationally broadened 
spectra resides in the metal lines and in central parts of 
hydrogen and helium lines. It is important to get rid of 
the continuum and of the broad wings of Stark-broadened 
lines when calculating the figure of merit, otherwise these 



components dominate the determination of usim, whereas 
in practice they carry no information on usim. This prob- 
lem can be solved by high pass filtering both spectra before 
comparing them. The filtering is done with a high-pass fre- 
quency filter in Fourier-space. The details of the filtering 
process are given in Appendix B. 

Figures 4 and 5 show an example of the fitting process 
for a star for which the best value of usim is 118 kms -1 . 
It is HD 121790, a B2IV-V, star and the template is 
HD 96706 with spectral type B2V. Figure 4 shows a 
direct comparison of the broadened template and the 
target spectrum. Figure 5 shows a comparison of the 
Fourier filtered spectra. All spatial wavelengths longer 
than 358 kms -1 were filtered out in this case. Note that 
in Fig. 4 there are several mismatches between the two 
spectra in the continuum as well as in the wings of the 
hydrogen lines. These mismatches bias the determination 
of usinz if no filtering is applied. The comparison between 
the two filtered spectra shows convincingly that the broad- 
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ened template spectrum fits the target spectrum. The in- 
terstellar Call A3934 line is located in the central part of 
order 145. This line is left out of the comparison between 
the spectra. 

As mentioned in the beginning of this section, the tem- 
plate broadening technique fails at high values of vsini. 
The reason is that the classical model for rotational broad- 
ening does not take the effects of gravity darkening into ac- 
count. For stars rotating at a large fraction of their break- 
up velocity the surface gravity decreases in the equatorial 
regions. This leads to obscuration of the parts of the ob- 
served line profile that come from the equatorial regions. 
The result is that rapidly rotating stars which are observed 
equator-on have line profiles that resemble those of more 
slowly rotating stars. Furthermore, because of the strong 
lines that are used we expect the effects of limb-darkening 
to be more severe. We investigated both these effects on 
the template broadening method by using synthetic data 
generated from the models by CTC. These investigations 



Fig. 5. The same as Fig. 4 but now 
for the high-pass filtered spectra. The 
order numbers and central wavelengths 
are the same as in Fig. 4. The distance 
between the tick-marks on the vertical 
scale is 0.01. The high-pass filter takes 
out all spatial wavelengths longer than 
358 kms -1 (see also appendix B). No- 
tice the very good fit of the two spectra. 

are described in Appendix C. The conclusions are that 
the choice of the limb-darkening law is not important for 
template broadening, but gravity-darkening is important 
and one should be careful when using template broadening 
results as soon as vsini is larger than ~ 120-150 kms~ . 

3.3. Determining vsmi for Rapid Rotators by Model 
Fitting 

For the stars with the highest values of usinz one is only 
left with the alternative of comparing model line profiles 
to the data. The only line in our spectrum for which rotat- 
ing stellar models have been calculated is the Hei A4026 
line. The models are given by CTC for a set of values of 
w, the fraction of the critical angular velocity at which 
the star rotates, and i, the inclination angle of the ro- 
tation axis with respect to the line of sight. The critical 
angular velocity of the star corresponds to the equatorial 
rotational velocity at break-up. The models were calcu- 



A.G.A. Brown & W. Verschueren: Rotational Velocities of Early- Type Stars in Sco OB2 



9 




4000 4010 4020 4030 4040 4000 4010 4020 4030 4040 

A (A, logarithmically spaced) 



lated for main sequence stars and for stars of luminosity 
class III for the spectral types Bl, B3, B5, B7, and B9. 
The CTC models are based on the Roche model described 
in Collins (1963, see also Sect. 4.2). The basic procedure 
we follow is to compare the whole line profile to the data 
for each (w,i) combination and determine vs'mi from the 
best fitting models. We stress here that we cannot derive 
w and i separately with our data. Usually a number of 
(w, i) combinations are consistent with the data and one 
can only derive vsini with a corresponding range of un- 
certainty. 

From the set of models given by CTC we constructed 
a grid of models for each spectral type. In this grid w 
varies from to 1.0 in steps of 0.1 and i varies from 0° to 
90° in steps of 10°. The models at grid points not given 
by CTC were constructed by linear interpolation between 
CTC models. We verified that using higher order inter- 
polation methods does not improve the model grid. We 
also constructed models for B4, B6 and B8 stars by lin- 
ear interpolation between spectral types given by CTC. A 
similar procedure is not possible for B2 stars because they 
have the strongest He I A4026 lines; and here interpolation 
between Bl and B3 models is incorrect. For B2 stars B3 
models were used and we discuss the implications below. 

The model profiles are given as a function of wave- 
length so we used spectra rebinned into In A-space. The 
He I A4026 profiles were transformed to In A-space (spac- 
ing in In A: 2.8 x 10~ 5 ). The data and the models were 
compared for each (w, i) pair and the pair for which the re- 



Fig. 6. Three examples of the model 
fitting process for deriving vsini. Fig- 
ure (a) shows a direct comparison of the 
best fitting model (thick line; w = 0.6, 
i = 50°, vsini — 130 kms -1 ) and the 
data for HD 108250. Figure (b) shows 
the same as (a) after high-pass filtering. 
Figures (c) and (d) show the same as (a) 
and (b) for HD 134481. The best fitting 
model in this case is: w = 1.0, i = 30°, 
vsini = 174 kms -1 . Figures (e) and (f) 
show the same as (a) and (b) and illus- 
trate the effect of using a B3 model to 
fit the data of a B2 star. The He I A4026 
line is weaker for B3 stars than for B2 
stars and the effect is that vsini is un- 
derestimated. The best fitting model in 
this case is: w — 0.5, i = 60°, vsini = 
124 kms -1 . The dashed lines indicate 
the continuum level of the models and 
illustrate the need for high-pass filter- 
ing. 

duced x 2 is a t a minimum was considered to correspond to 
the best fitting model. Possible radial velocity shifts were 
taken into account by automatically shifting the model 
with respect to the data for each (w, i) pair and searching 
for the minimum \ 2 . 

In order to estimate the uncertainties in the derived 
value of vsini the ten models with the lowest x 2 were 
considered. Of these, a model was considered equivalent 
to the best fitting model if the associated x 2 did not dif- 
fer by more than a factor of 2 from the minimum and if 
the emergent flux at all points in the model did not dif- 
fer by more than 0.01 for Hci A4026 lines with a central 
depth of 0.9 or less. The maximum tolerated difference in 
fluxes was 0.005 for Hci A4026 lines with a central depth 
of more than 0.9. The vsini values corresponding to the 
best fitting models were subsequently averaged with the 
corresponding x 2 values as weights. The rms spread in 
vsini values was taken as the error bar. In some cases one 
finds only one best fitting model and then the adjacent 
points located between the grid points were considered in 
order to assign an error bar to vsini. 

As was the case with template broadening, there may 
be mismatches between model and data due to the contin- 
uum. We turned to high pass filtering again to get around 
the continuum terms. Because one wants to compare the 
whole line profile to the data, and because most of the 
target stars have a high value of vsini, only the most im- 
portant continuum terms should be filtered out. In or- 
der to safely examine stars with values for vsini up to 



10 



A.G.A. Brown & W. Verschueren: Rotational Velocities of Early- Type Stars in Sco OB2 



500 kms -1 , we chose to filter away all spatial wavelengths 
larger than ~ 1000 krns" 1 . None of the stars in our sam- 
ple have vsini values this high. 

Figure 6 shows three examples of the model fitting pro- 
cess. In all cases the comparison between data and best 
fitting model is shown both for the high-pass filtered data 
and for the un-filtered data. In Figs. 6a and 6b the data 
for HD 108250 (B4V) are shown. The best fitting model 
is indicated by the thick line. Note that in the direct com- 
parison it appears as if the Hei A4026 model line is not 
deep enough and that the wings of the line profile do not 
match the data either. This is due to the continuum in the 
data which has not yet been properly rectified. This can 
be seen from Fig. 6b; the Fourier filtered model and data 
now match well in both the line core and in the wings. 
Figures 6c and 6d show the same for HD 134481 (B9V). 
In this case the continuum close to the Hei A4026 line just 
happens to have the right value. 

The strength of the He I A4026 line reaches a maximum 
around spectral type B2-B3 and the B3 models are not 
always strong enough compared to the data. This is often 
the case when one considers B2 stars. Figures 6e and 6f 
show the effects of using B3 models to fit the data for a 
B2 star (here HD 129116 B2.5V). The Hei A4026 line is 
weaker in the B3 model than in the B2 star. Because both 
the broader wings and the deeper line core of the B2 star 
have to be fitted, vsini is underestimated. 

All results were inspected by eye to judge whether the 
best fitting model was actually correct. We encountered 
stars for which the models corresponding to their spectral 
types did not fit the data. In those cases other models 
were tried and the spectral type for which the best fit 
was obtained was used to derive vsini. For the majority 
of these cases the strength of He I A4026 and the ratio of 
Hei A4026 to Hei A4009 in our spectra indeed were not 
consistent with the spectral type of the star. 

Finally, we remark that the vsini values are based on 
the rotational velocity of the star at the theoretically cal- 
culated break-up velocity, and the value of w and i. This 
means that the later interpretation of vsini distributions 
in terms of break-up velocities will be complicated (see 
Sect. 4). 

3.4- Results 

The results of the vsini determinations are shown in Ta- 
bles 2a-2d. The tables contain the following data. Col- 
umn (1) gives the HD number of the star and column (2) 
gives the number of the star in the HIPPARCOS Input 
Catalogue (Turon et al.1992). In column (3) the spectral 
type is given as listed in the HIPPARCOS Input Cata- 
logue. The derived value of the projected rotational veloc- 
ity and its associated error are listed in columns (4) and 
(5). For stars for which we have vsini values from different 
runs the error-weighted average is listed. For the following 
sharp-lined stars the vsini values from two or more runs 



are not consistent within the error bars: HD 109668, HD 
110879, HD 120307, HD 122980, HD 132955, HD 133955, 
HD 145502, HD 147165, HD 158408, HD 172910. These 
stars all turn out to be spectroscopic binaries or radial ve- 
locity variables (see below). CTC provide a value for the 
break-up velocity («b r ) for each spectral type for which 
model profiles were calculated. Columns (6) and (7) list 
this break-up velocity and the ratio vsini/v^. In cases 
where the CTC model profiles were used the value of «b r 
corresponds to the spectral type of the best fitting model 
and not necessarily to the listed spectral type (see also 
the notes to Table 2). For stars later than AO and for su- 
pcrgiants we determined the break-up velocities based on 
stellar parameters given by Straizys & Kuriliene (1981) 
and Harmanec (1988). For HD 149757 and HD 151804 we 
took the mass and radius given by Lamers & Leitherer 
(1993). In column (8) we list which method was used to 
derive vsini. The method is indicated by its number and 
for some stars more than one method was used and a 
weighted average of vsini was then formed. For a num- 
ber of stars none of the methods for determining vsini 
works (these are very early-type stars or late-type stars 
with only blended lines, or spectroscopic binaries for which 
line-profile fitting does not work). For these stars vsini was 
estimated by visual inspection of the spectrum or by tak- 
ing a literature value; as indicated in the notes. In column 
(8) this is indicated by the number 4. Column (9) indi- 
cates whether or not a star is a member of Sco OB2; 'm' 
indicates that the star is a photometric member, as found 
by de Geus et al.(1989), and 'M' indicates that the star is 
a proper motion member, as found by Blaauw (1946) and 
Bertiau (1958). Finally, in column (10) we list whether 
the star has a constant radial velocity (CON), is a radial 
velocity variable (RV), a single- lined spectroscopic binary 
(SB1) or a double- lined spectroscopic binary (SB2). For 
stars without an entry in column (10) there is no good 
information on their binarity. 

The duplicity information was obtained from the 
Bright Star Catalogue (Hoffleit & Jaschek 1982; Hoffleit, 
Saladyga & Wlasuk 1983), Levato et al.(1987) and from 
preliminary findings by Verschueren et al.(1996) based on 
our data. A star is considered a radial velocity variable 
if it is listed as SB or SB1 in the Bright Star Catalogue 
or if Levato et al.or Verschueren et al.list it as a radial 
velocity variable. A star is considered a single-lined spec- 
troscopic binary if the Bright Star Catalogue lists it as 
SBO or SBIO or if Levato et al. determined an orbit and 
list it as SB1. A star is considered a double- lined spec- 
troscopic binary if any of these references list it as SB2 
or SB20. Note that the duplicity information pertains to 
duplicity that is apparent in the spectrum of the observed 
object; either through variability of the radial velocity or 
because of the presence of double lines. Thus, although 
some stars are part of a (wide) binary or multiple system, 
if there are no signs of duplicity in their spectrum they 
are listed as single stars. Finally, an asterisk after column 
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Table 3. Projected rotational velocities for early-type 
stars not related to Sco OB2. 
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(10) indicates that there is a note about that star in Ap- 
pendix D. The notes contain information on the derivation 
of vsini. They also provide information about non-radial 
pulsations in the stars. These pulsations can have a size- 
able effect on the line profiles and thus on the derived 
value of vsini. The information on the non-radially pulsat- 
ing stars in our sample (including the (3 Cephei variables) 
was kindly provided by P. Groot. 

For the star HD 140008 it was possible to measure 
vsini for both components in the spectrum, and to within 
the error bars they both have the same vsmi of 11 kms -1 . 
Both components were taken into account in the subse- 
quent analysis of the data. For method 2 the following 
stars were used as templates: HD 96706, HD 120709, HD 
126341, HD 129056, HD 132955, HD 134687, HD 149438. 
In Table 3 we give the projected rotational velocities of 
the additional set of early-type stars that was observed 
with ECHELEC. Of these stars the following have been 
used as templates for method 2: HD 28114, HD 35039, HD 
36591, HD 45572, HD 72660, HD 179761. 

Finally, the stars HD 140008, HD 143118 and HD 
144294 are listed as members of UCL by de Geus et al.(1989) 
but in order to be consistent with the subgroup bound- 
aries as given by de Zeeuw et al.(1994) we list these stars 
as members of US. 

3. 5. Consistency of the Three Methods for Determining 
vsmi 

Figure 7 shows a comparison of the vsmi values derived 
using the three methods described above. In Fig. 7a the 
results from line- width measurements and template broad- 
ening arc compared. All the values are consistent within 
the error bars. However, the values from line-width mea- 
surements are systematically lower, although generally 
this effect is ^3 kms -1 . This is an effect of the lowering 
of the apparent continuum due to rotational broadening. 
A lower continuum means that the measured line-widths 
are too small, which leads to underestimates of vsmi. 

Figure 7b shows the comparison of vsmi values de- 
rived from template broadening and fitting the He I A4026 
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Fig. 7. Comparison of vsini values obtained from measure- 
ments of line-widths, template broadening and fitting models 
to the He I A4026 line. Figure (a) shows the comparison be- 
tween the line-width method and template broadening. Figure 
(b) shows a comparison of the results from template broad- 
ening and model fitting and (c) shows the same for stars of 
spectral type earlier than B5. The B2 stars are indicated as 
black dots. Figure (d) shows the same as (c) but now the tem- 
plate broadening values are derived from order 142 (containing 
the He I A4026 line) only. See the text for discussion. 

line with models. Template broadening was applied for all 
spectral types over a large range of vsini. Thus, bearing 
in mind the discussion in Sect. 3.2 and Appendix C, we 
ilxpect the two methods to yield inconsistent results for 
vsini 150 kms -1 . This can be clearly seen in Fig. 7b. 
For 17 out of the 52 stars in this figure the vsini values 
from each method are inconsistent if one takes the error 
bars into account. 

In order to understand Fig. 7b better we break the plot 
down according to spectral type. For spectral types earlier 
than B5 most stars have vsini/vbr r$0.5. The opposite is 
the case for later spectral types. So for later spectral types 
w > 0.7 (see Sect. 4 for a discussion on the relation be- 
tween vsini/vbr and w). That is precisely the range where 
template broadening fails. Figure 7c shows the compar- 
ison of the two methods if one leaves out spectral type 
B5 or later. The scatter at the high vsini end is much 
reduced but overall it is still significant. We have also in- 
dicated the B2 stars separately and one can see that most 
of them have model fitting values that are lower. Wc noted 
in the previous section that this is to be expected if the B3 
He I A4026 model lines are not strong enough compared to 
the data. 
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Figure 7d shows the comparison if one uses only the 
results for template broadening from order 142 (contain- 
ing the He I A4026 line). Most of the B2 stars are displaced 
downwards, which is consistent with the remarks about B3 
models made above. The other stars now lie closer to the 
line of consistent values and overall the scatter is some- 
what reduced. The results in this figure suggest that the 
Balmer lines in other ECHELEC orders lead to underes- 
timates of vsini for template broadening. This is indeed 
confirmed if one compares the template broadening results 
of all other ECHELEC orders to those from order 142. 

Another cause for inconsistencies between template 
broadening and model fitting comes from the fact that 
one has to make assumptions about the stellar parame- 
ters in model fitting. One needs to know the radius and 
mass of the stars in order to calculate vt, r and hence vsini. 
This problem does not arise when one uses the classical 
model for rotational broadening. For example: CTC give a 
mass of 11 M Q and a radius of 5.8 R Q for a B1V star. How- 
ever, according to the calibrations by Straizys & Kuriliene 
(1981) a B1V star has a mass of 12.9 M and a radius of 
5.9 R , while Schaller et al.(1992) use 9 M Q and 3.7 R . 
These parameters lead to break- up- velocities of 491, 527 
and 556 kins" 1 , respectively. So a variation of more than 
10% in vsini is possible due to uncertainties in stellar pa- 
rameters. 

We conclude that the three methods for obtaining 
vsini are on the same velocity scale. However, caution 
should be exercised when comparing the results from tem- 
plate broadening to those of model fitting. The average 
difference between line- width and template broadening re- 
sults is 7%, and for template broadening and model fitting 
this difference is 10%. 
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Fig. 8. Comparison of the vsini values derived by SCBWP and 
the values derived in this paper. 



3.6. A Comparison of vsini Values for Standard Stars 

We now turn to the comparison between our results 
and the values for vsini listed for the standard stars in 
the system presented by Slettebak et al.(1975, hereafter 
SCBWP). They determined projected rotational velocities 
for stars with a range of spectral types and vsini. These 
measurements define a vsini velocity scale to which other 
measurements can be compared. 

Figure 8 shows the comparison between vsini de- 
rived by us and the values listed in SCBWP (1975). We 
have 51 stars in our sample in common with SCBWP. 
The error bars on the values of SCBWP are assigned 
according to the accuracy given in their paper (10% 
for vsini ^200 kms -1 and 15% for vsini ^200 kms" 1 ), 
where for values of vsini less than 50 kms -1 an error bar 
of 5 kms -1 is assigned. Within the error bars almost all 
vsini values are consistent, but below about 60 kms -1 
most of our values are systematically offset from those 
by SCBWP to higher vsini. The values of SCBWP are 
based on model calculations presented by Collins & Son- 
neborn (1977). They compared the equivalent width ratio 
of He I A4471 to Mgn A4481 to the model calculations to 
derive vsini. This procedure is not easy to compare to 
what we have done. 

However, if we use the line-widths for Mgn A4481 as 
measured by SCBWP and calculate vsini according to 
the modern models of CTC the discrepancies largely dis- 
appear and are no longer systematic. The new vsini val- 
ues were derived directly from the Mgn A4481 line- widths 
according to the procedure described in Collins & Truax 
(1995) after correction for the intrinsic and instrumen- 
tal broadening of the line. Among the sharp-lined stars 
there are also two giants for which SCBWP used main- 
sequence models to derive vsini, which leads to underes- 
timated values. Therefore, we conclude that our derived 
values of vsini are on the same velocity scale as defined 
by SCBWP. 

4. The Distribution of Rotational Velocities in 
Sco OB2 

In this section we discuss the distribution of projected 
and true rotational velocities in Sco OB2. We concentrate 
on the known members of the association. Of the 156 
stars we observed 88 are known members of the associ- 
ation (counting the double-lined spectroscopic binary HD 
140008 as two objects). In order to increase the sample size 
we added vsini values from literature for an additional 48 
members (counting HD 145501 twice). The membership 
of these stars was taken from de Geus et al.(1989) and 
Bcrtiau (1958). The vsini values were taken from Slette- 
bak (1968) and Uesugi & Fukuda (1981). Their values are 
on the old velocity scale defined by the studies of Slette- 
bak (1954, 1955, 1956, 1966a, 1966b, 1968) and Slettebak 
& Howard (1955), whereas our vsini values are on the 
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velocity scale defined by SCBWP. In order to take the 
differences into account we corrected the literature values 
downwards by 15% and 5% for B-type stars and A-type 
stars, respectively, as recommended by SCBWP. In Table 
4 we give the resulting values of vsini for the literature 
sample. The columns have the same meaning as in Table 
2. The stars HD 120908 and HD 138690 are members of 
UCL, all others arc members of US. Note that HD 147889 
is a doubtful member according to Bertiau (1958) and de 
Geus et al.(1989) do not list this star as a member. How- 
ever, including or excluding this star from the analysis 
does not influence the results. 

Table 5 shows the distribution among spectral type 
for the members of the subgroups of Sco OB2. The num- 
bers in parentheses indicate the stars for which vsini was 
obtained from literature. Note that the sample observed 
by us is very biased towards spectral types B0-B3 (see 
Sect. 2). If one assumes the initial mass function given in 
de Geus (1992) for Sco OB2, the numbers indicate that 
our observations are far from complete. Note also that of 
the 48 vsini values from literature 46 are for stars in US 
and of these 33 are of spectral type B7-B9. All this means 
that caution should be exercised when comparing distri- 
butions of vsini for the various subgroups. 

4.1. The vsini Distributions in the Sco OB2 Subgroups 
and the Binary Population 

Figures 9 and 10 show the distributions of vsini and 
vsini/vb r for Sco OB2 and each of the subgroups for which 
membership studies have been done. In both figures we in- 
dicate the distributions of single stars and binaries (spec- 
troscopic binaries and radial velocity variables) separately, 
with solid and dotted lines respectively. Both figures sug- 
gest that amongst the binaries there is a higher proportion 
of slow rotators. We tested this by performing a robust- 
rank-order test on the members of Sco OB2 that we have 
observed and for which we have good information on du- 
plicity (see Table 2). This is a non-parametric statistical 
test that is designed to test for differences in the median of 
two samples (see e.g., Siegel & Castellan 1988). Nothing is 
assumed a priori about the underlying distributions. For 
the vsini distributions this test confirms at a confidence 
level of more than 99.99% that the single stars have a 
higher median vsini. For the v sini/vbr distributions the 
same result is obtained. This implies that the bulk of the 
single stars has a higher projected rotational velocity. We 
find the same result if we use all the members of Sco OB2 
(assuming that stars for which we have no duplicity infor- 
mation are single). 

The finding that binaries show lower rotational veloci- 
ties can be interpreted as a consequence of the deposition 
of part of the total angular momentum in the orbit of the 
system at formation or as a consequence of tidal braking in 
very close binaries. Levato ct al.(1987) have shown, using 
the vsini data of Slettebak (1968), that there is a corre- 



lation between the frequency of duplicity and the average 
vsini for the subgroups US and UCL. This correlation is 
in the sense that stars in UCL rotate slower on average 
and that the frequency of duplicity is higher in that sub- 
group. Based on the duplicity data in Tables 2 and 4 we 
also find a higher proportion of binaries in UCL than in 
US. However the distributions of rotational velocities are 
not significantly different for these two subgroups. 

If one looks at the data in more detail it turns out 
that most of the single stars at high vsini values are stars 
added from the literature in the spectral type range B7- 
B9, whereas our sample contains mainly B0-B6 stars (see 
Sect. 2). Most of the literature data (46 out of 48 stars) was 
added to US and indeed the difference between the single 
and binary stars is most pronounced for this subgroup (see 
Figs. 9 and 10). In US the extra literature data consist of 
28 stars for which vsini > 150 kms -1 . Of these stars 15 
are single stars in the B7-B9 range. The sample observed 
by us contains, in contrast, only one single B7-B9 star 
in the same vsini range. To check whether the stars in 
the B7-B9 range rotate faster we divided the sample in 
B0-B6 and B7-B9 single stars. Each of these samples was 
then divided into stars for which vsini > 150 kms" 1 and 
vsini < 150 kms" 1 . Subsequently, a 2 x 2 contingency 
table analysis was done on the data with the Fisher exact 
probability test (see e.g., Siegel & Castellan 1988). The 
test confirms at a confidence level of more than 99.99% 
that the B7-B9 single star sample contains more objects 
with vsini > 150 kms" 1 . If one does the same test for 
vsini/vbr with the division of the samples at 0.4 then the 
same result is found at the 99.98% confidence level. The 
value 0.4 follows from the fact the all B7-B9 stars with 
vsini > 150 kms" 1 have vsini/vb r > 0.4. Thus, it could 
be that late B-stars exhibit intrinsically faster rotation. 
However, this finding is based on inhomogeneous vsini 
data. To really interpret the vsini data, a more complete 
coverage of the spectral types is needed in all subgroups of 
Sco OB2 and a complete census of the binary population 
as a function of spectral type is also required. 

Finally, we remark that there are no significant differ- 
ences between the subgroups of Sco OB2 for the median 
value and the overall distribution of vsini or vsini/vbr- 

4-2. Tests of Various Hypotheses about the Distribution 
of Intrinsic Rotational Velocities 

We now test various hypotheses about the distribution of 
intrinsic rotational velocities against the data. We wish to 
investigate whether any constraints can be placed on the 
distributions of true rotational velocities and/or the ori- 
entation of rotation axes. This will lead to constraints on 
the origin of angular momentum in stars and on the star 
formation process. So far the only models for vsini distri- 
butions that provide detailed predictions are the ones that 
assume a random selection process. The magnitude of the 
rotational velocity and the inclination angle of the rota- 
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Table 4. Projected rotational velocities for the literature 
sample. 
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^^^lion axis are selected separately and randomly according 

to a certain prescription. These models do not take into 

account the physically more likely possibility that a star 
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Table 5. Distribution of spectral types for the members 
of Sco OB2. 



Subgroup 





B0-B3 


B4-B6 


B7-B9 
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US 
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25 (5) 


3(7) 


3 (33) 


1(1) 


UCL 





18 (1) 


6(1) 


5 





LCC 





18 


6 


2 






is formed with a certain angular momentum. The role of 
stellar mass and the star formation process is also ignored. 

We assume in all cases that the rotation axes are ori- 
ented randomly in space. This means that the distribution 
of inclination angles i is given by: 



f(i) di — sinidi . 



(1) 



In what follows we define v c to be the true equatorial rota- 
tional velocity of the star, rj — w c /wbr and y = vsini/vbr- 
Again, «br is the true equatorial velocity of the star at 
break-up, and w is the fraction of the corresponding crit- 
ical angular velocity Note that in case the Roche model 
is assumed for the rotating star, vsini is not equal to 
«br x w x sini and that it makes a difference whether 
rj or w is chosen randomly to generate the distribution of 
y. In the Roche model the star rotates uniformly and its 
surface follows an equipotential surface. The potential is 
generated by a point mass at the centre of the star. The 
equation describing the surface is: 



GM 



1 



> 2 i? 2 sin 2 9 



GM 
R 



(2) 



$ is the total potential consisting of a gravitational and 
a centrifugal term. The mass of the star is M and the 
distance from the surface to the centre is given by R. The 
polar radius is i? p , which remains constant, and 9 is the 
latitude on the surface of the star as measured from the 
pole. The critical angular velocity follows by setting g = 
-V* = and is given by: w 2 = 8GM/27i?3. With w = 
w/u c and setting x = R/Rp Eq. (2) becomes: 



1 4 2 2 • 2 

1 w x sin i 

x 27 



1 . 



(3) 



Now, from the definitions of u> c , x and w it follows that 
rj = |u>x, where x is to be evaluated at the equator. This 
value of x follows by setting 9 — tt/2 in Eq. (3) and solving 
this equation for x for a given w. The equation can be 
solved analytically and it follows that: 



i] = —2 cos 



arccos w 



2tt 



(4) 



It is this transformation between to and r\ that one has to 
take into account if the distribution of y is generated by 
picking w randomly. 

We tested the following models for the distribution of 
intrinsic rotational velocities of the stars: 



1. The value of w is distributed homogeneously on [0,1]. 

2. The value of rj is distributed homogeneously on [0, 1]. 

3. The distribution of v c is given by a Maxwellian. 

4. The distribution of w is given by a truncated Maxwellian.| 

All the hypotheses about the distribution of v e are vari- 
ations on the assumption of a completely random distri- 
bution of rotational velocities. So the question is whether 
one can exclude any form of total randomness in the v e 
or w distribution. If all these models are excluded by the 
data then either one has to relax the assumption of ran- 
dom orientation of rotation axes or a non-random process 
for generating the distribution of rotational velocities has 
to be assumed. We find no correlation in our data between 
vsini and Galactic longitude, which is consistent with the 
assumption of random orientation of rotation axes. 

For the four models above one can derive an expression 
for the distribution of y by using the formula given by 
Chandrasekhar & Munch (1950): 



0(vsini) = vsini 



f(Ve) 



v c ^Vc 2 - (vsini)'- 



: dv e , 



(5) 



where 0(vsini) is the distribution function of vsini and 
/(v e ) that of v c . If there is a break-up velocity, then Eq. (5) 
becomes (in terms of y and rj): 



4>{y) = V 



y v\nf 



y 



:dr), 



(6) 



The form of 4>(y) for model 2 is derived in Chandrasekhar 
& Munch (1950, their Eq. (26)). To derive 4>(y) for model 
1 one should use 



with: 



. dw . 



(7) 



(8) 



Dcutsch (1970) showed, using arguments from statistical 
mechanics, that if u> is distributed isotropically with mu- 
tually independent Cartesian components the distribution 
of lo is given by a Maxwellian. The form of the distribution 
function is: 



f(z) = ^z\~ 



(9) 



where z = ju; and j is a parameter with the dimensions 
lo^ 1 . It follows that the distribution of vsinz is then given 
by: 



<j>{u) = 2ue~ u , 



(10) 



where u is jvshii and 1/j is the root mean square value 
of vsini which is the same as the modal value of v c for a 
Maxwellian (see Deutsch 1970). This is the distribution of 
vsini for model 3. Note that it is assumed that there is no 
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Fig. 11. (a) The distribution functions 
of y = vsini /v hr for models 1, 2 and 4. 
(b) The distribution function of usim 
for model 3. The models are indicated 
by their corresponding numbers. 



break-up velocity in this case. To derive the distribution 
for model 4 one should take into account that u < uu c . We 
chose to model this simply by truncating the Maxwellian 
at lv c and normalizing the resulting distribution. In terms 
of r) one can find <fr(y) by using Eqs. (4), (6) and (7). The 
distribution of w is now given by: 



w 2 e 



erf(l)">- 



(11) 



where erf is the error function. 

The distribution functions for y for models 1, 2 and 4 
are listed in Table 6 (model 3 being given by Eq. (10)). 
For model 4 the integral in Eq. (6) has to be evaluated 
numerically. All the 4>{y) are shown in Fig. 11. Note that 
models 3 and 4 show a peak in the distribution at inter- 
mediate values of vsini and vsmi/vbr whereas models 1 
and 2 predict a large number of slow rotators. 

Table 6. Distribution functions for y = vsmi/vb r . 



Model 



~2 y V* y2 + 2 arccos 



arccos y 



4 ^(l-^ 2 ) ^ 



erf(l)- w e 1 J v ^rf - y 



■. drj 



w = cos ( arccos (-^-)) 



The models were compared to the data by calculating 
the Kolmogorov- Smirnoff statistic D and the probability 
P that D exceeds the observed value under the hypothesis 
that the data were drawn from the model distribution (see 



e.g., Press et al.1992). A model was rejected if P < 0.05. 
As argued in the previous section, the B7-B9 stars may 
be a group of stars that rotate intrinsically faster than 
B0-B6 stars. Our data consist primarily of B0-B6 stars 
and adding literature data on later types only to US may 
introduce artifacts in the vsini distribution that are the 
result of our sample selection. Thus we decided to compare 
the model distributions only to the B0-B6 stars. 

We did the KS-test for Sco OB2 as a whole and for each 
subgroup. The data were also divided in single stars and 
binaries. The binaries include the spectroscopic binaries 
as well as the radial velocity variables. The comparison of 
the models to the data was focused on the members of Sco 
OB2. We are interested here in the distribution of rota- 
tional velocities in a group of physically related stars. The 
results are shown in Table 7. In this table a "+" indicates 
that the model was not rejected and a "— " indicates that 
the model was rejected. 

The number of measurements for each subgroup is 
rather small (40 stars in US, 26 in UCL and 24 in LCC), 
which makes it difficult to distinguish between models. 
The data for the single stars are compatible with all mod- 
— els. This is probably the result of the limited amount of 
ydata. If we take only the data for US into account and 
include the B7-B9 stars we find that only model 4 is con- 
sistent with the data. It is not warranted, however, to 
conclude that the formation of single stars leads to an 
isotropic distribution of angular momentum vectors. We 
have shown in the previous section that the B7-B9 stars 
rotate faster than B0-B6 stars. So there is an effect with 
stellar mass that model 4 does not take into account. 

For the binary stars it seems that only model 1 is con- 
sistent with the data. What does it mean if we assume that 
this model actually describes the data? Do components of 
—binary stars acquire a random fraction of their critical 
angular velocity or is model 1 the only model that pre- 
dicts enough slow rotators? Remember that binary com- 
ponents show intrinsically slower rotation than the single 
stars. This may be due to the redistribution of angular 
momentum into orbital angular momentum. 
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Table 7. Results of comparing model distributions to the 
data (B0-B6 members). A "-" indicates that the model is 
rejected. A "+" indicates that the model is not rejected. 
The number of stars in each data set is indicated in paren- 
theses. 
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If one aims at putting real constraints on the star 
formation process by looking at the distribution of ro- 
tational velocities, then models arc needed that take into 
account the details of this process. How does the com- 
bination of stellar mass, angular momentum and interior 
structure translate into an observed rotational velocity? 
So far no detailed predictions of this kind have been made. 
Mouschovias (1983) does give an estimate of the equato- 
rial velocity of the star as it arrives on the zero-age main- 
sequence. His prediction follows from a consideration of 
magnetic braking in the early stages of star formation. As 
soon as the protostellar condensation is dense enough to 
decouple from the surrounding magnetic field its angular 
momentum is basically determined. It then follows that: 



oc p 2 / 3 RM c ■ 



(12) 



where R and p are the mean density and radius of the pro- 
tostar and Mq is the mass of the cloud-core from which the 
protostar formed. The factor p 2 / 3 R increases going from O 
to F0 stars. If one assumes that the mass of the cloud-core 
scales with the mass of the star that forms from it, this 
equation predicts a maximum in the rotational velocities 
somewhere between O and F stars. This is consistent with 
the observed distribution of velocities for field stars (e.g., 
Wolff et al. 1982), and with our observed distribution of 
rotational velocities as a function of spectral type. 

Another model for the origin of angular momentum in 
stars that takes the star formation process into account is 
given by Wolff et al.(1982). They find a large number of 
slow rotators in their data that cannot be explained by the 
presence of binaries. Therefore they propose that all pro- 
tostars start out in a state of low angular momentum. Sub- 
sequently, a number of them acquire higher angular mo- 



menta through gravitational interactions with other pro- 
tostars or protostellar condensations in a forming cluster. 
This model then predicts a concentration of rapid rotators 
in the centres of clusters and associations. We cannot see 
this clearly in our data and in our case the presence of 
binaries can account for the excess of slow rotators with 
respect to a Maxwellian. 

To go further more detailed predictions of the distri- 
bution of rotational velocities are needed that take the 
star formation process into account. Observationally, more 
data are required (i>sim as well as duplicity information) 
for the members of Sco OB2 of spectral types later than 
B6. This will facilitate a study of usim as a function of 
mass and of the relation to duplicity. 



5. Rectifying the ■usini Distribution 

In the previous section we compared model v e and v e /vhr 
distributions to the observed wsini distribution in Sco 
OB2. In this section we present the distributions of v c and 
v e /t>br that follow directly from the rectification of the ob- 
served vsini distribution under the (sole) assumption of 
randomly oriented rotation axes. Two techniques are used 
to do so; Lucy-iteration (Lucy, 1974) and the technique 
presented by Bernacca (1970). 
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Fig. 12. Results of applying the Lucy iterative technique for 
the rectification of usini distributions. The number of itera- 
tions used was 8, 9, 11 and 9 for LCC, UCL, US and the ECH- 
ELEC data in US, respectively. The data histograms are shown 
for comparison. The solid lines indicate (p(vsmi), the dashed 
lines are for tp(v e ) and the short-dashed lines are for y>(usini). 
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For both techniques it is necessary to have a good esti- 
mate of the observed probability density function (PDF). 
We used a so-called kernel estimator. The PDF is repre- 
sented as a sum of kernel functions, each corresponding 
to a single data point. The kernel function is normalized 
to unity and it is characterized by a width h. The width 
depends on the observed distribution of the data. A rule 
of thumb for calculating h is given in Vio et al.(1994). The 
kernel estimator provides a smooth estimate of the PDF. 
We used Gaussian kernels, in which case h corresponds to 
the standard deviation. 

For both rectification techniques the interval for the 
kernel-estimate of the PDF was restricted to vsmi > 
kms -1 . For the Lucy-iteration the upper limit of the in- 
terval was set by increasing the data interval in small steps 
until the cumulative distribution increased by no more 
than 0.01%. In the technique given by Bernacca (1970) 
one assumes that there is an upper limit to usini. We 
used the break-up velocity of the stars as a function of 
spectral type and rectified the vsvai/vbr distribution. So 
the PDF was restricted between and 1. In both cases the 
kernel-estimator formally has tails outside the intervals 
mentioned above, and we decided to simply re-normalize 
the PDF on the restricted intervals. 

For the iterative technique given by Lucy (1974) we 
used the optimum stopping criterion described in Lucy 
(1994). In his paper two quantities are defined: 



N 



H = ^2 & 111 <Pi 



and 



i=i 



N 



(13) 



(14) 



where H is A^ 1 x the log-likelihood of the data vector 
(p (this vector being the estimate of the PDF described 
above), with ip the data vector that one derives from the 
rectified distribution ip. S is a function that decreases as 
the solution ip increases in complexity (Note that we took 
Xj constant in Eq. (6) of Lucy 1994). The solution ip will 
describe a trajectory in the (H, S) plane during iteration 
and one can define the curvature of this trajectory as: 



k = 



SH-HS 



(S 2 - H 2 ) 



(15) 



It turns out that the iterations should be stopped when k 
passes through a minimum. 

Figure 12 shows the results of rectifying the observed 
vsmi distributions for the members of each subgroup of 
Sco OB2. For US we also show the result of rectification if 
we use only our own data. In all figures the data histogram 
is shown. The solid line is the estimate of the observed 
PDF (p(vsmi), the dashed lined is the rectified distribution 



ip(v e ) and the short-dashed line is the PDF ip(vsini) that 
follows from ip. 

Figure 12 suggests that the peak in the v c distribu- 
tion lies at increasingly higher values going from UCL to 
LCC to US. This result is the same even if only our own 
data for US are used. Whether this effect is due to a de- 
crease in binary frequency per subgroup is hard to tell. 
The data are incomplete for the B7-B9 stars for UCL and 
LCC. If one takes only B0-B6 stars into account and only 
our sample the binary frequency increases from LCC to 
UCL to US, opposite to the expected trend. However if 
one takes all the data into account then the binary fre- 
quency increases from LCC to US to UCL. We stress here 
again that the differences between the subgroups are not 
statistically significant. An improved sample is needed to 
settle this issue. 

It is interesting that the shift in the peak of the v e dis- 
tribution correlates with the ages of the subgroups. This 
suggests that the shift could also be due to evolutionary 
effects. Stars that evolve off the main sequence expand 
their outer layer and this may lead to a lowering of v c . 
This effect should be most pronounced in the giants and 
supergiants. The percentage of stars in these luminosity 
classes is, however, very similar for the three subgroups. 
If one also takes sub-giants into account then indeed UCL 
has the highest percentage of evolved stars, but the dif- 
ference with the other subgroups is small. Moreover UCL 
and US have similar percentages of evolved stars. Apart 
from the statistics, the magnitude of evolutionary effects 
depends much on whether or not stars rotate differentially. 

A drawback of the Lucy-rectification results presented 
above is that the existence of a break-up velocity is not 
taken into account. Yet, the wsinz values we derived for 
rapid rotators are based on the value of the break-up ve- 
locity for the star in question. The rectification technique 
presented by Bernacca (1970) does take the existence of 
an upper limit on v c into account. If there is a break-up 
velocity then one can define an angle 9 by: 



vsim — Wbr sint/. 



(16) 



It is clear that the inclination angle i must then satisfy 
9 < i < ir/2. So one has more information about the 
rotation axes of the stars. If one assumes that the rotation 
axes are randomly distributed in space with the restriction 
on i given above, the distribution of i, f(i), can be derived 
from the observed vsim distribution: 



/w = sm Vo n 



<Kv) 



r 



dy. 



The distribution of v c /vbr (rj) is then given by: 

y 2 <t>{y) 



i p 



V 7 ! 2 - y 2 V 1 - v' 



■dy. 



(17) 



(18) 



For a derivation see Bernacca (1970). We have translated 
his formulas to y and i]. 
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Thus, with the technique of Bernacca one can derive 
the distribution of i as well as rj and the distribution of 
■q is now given by an integral which can be evaluated nu- 
merically. The results for the three subgroups of Sco OB2 
are shown in Fig. 13. We used the same estimates of the 
observed PDF as for the Lucy-rectification as input for 
Eq. (18). The resulting distributions of rj (Fig. 13a) show 
the same trend as the v c distributions in Fig. 12. Namely, 
a shift of the peak towards higher values of rj going from 
UCL to LCC to US. Note that the up-turns in the distri- 
butions for LCC and UCL at 77 = 1 are probably not real 
but due to the singularity in the integral in Eq. (18). The 
distributions of i (Fig. 13b) are all very similar and essen- 
tially the difference from the distribution given in Eq. (1) 
is that there are less stars seen pole-on and more seen 
equator-on. 



6. Effects of Rotation in the Hertzsprung-Russell 
Diagram 

Rapidly rotating stars have deformed surface layers and 
their surface gravity varies from the pole to the equator. 
The atmospheric zones around the pole have close to nor- 
mal surface gravities and thus emit an almost normal stel- 
lar spectrum. The zones close to the equator, however, 
have lower surface gravities and also suffer from obscura- 
tion due to gravity darkening. The stellar spectrum emit- 
ted from these zones then corresponds to lower surface 
gravities as well as lower effective temperatures. This im- 
plies that one expects changes in the photometric colours 
of the stars with rotation. These changes are a function of 
the inclination angle as well as the rotation speed. 

These effects have been calculated by a number of 
authors for various photometric systems (see e.g., CTC 
and references in Tassoul, 1978). In general stars are 
displaced away from the main-sequence and the main- 
sequence for rotating stars always appears brighter than 
its non-rotating counterpart. Stars observed pole-on are 
displaced vertically towards higher luminosities and stars 



Fig. 13. (a) shows the distributions 
of 77 = Vc/vbr for the subgroups of 
Sco OB2 derived with the technique 
of Bernacca (1970). The solid line in- 
dicates the result for LCC, the short- 
dashed line that for UCL and the dashed 
line indicates the result for US. The up- 
turns in the distributions for LCC and 
UCL at rj = 1 are probably not real 
but due to the singularity in the integral 
in Eq. (18). (b) shows the correspond- 
ing distributions of i. In this figure the 
long-dashed line indicates the distribu- 
tion f(i) = sini. 

seen equator-on are displaced more along the main se- 
quence towards lower T c g ■ 

The effects of rotation on the Walraven photometric 
system have not been studied previously. In this section 
we investigate whether there are stars for which the Wal- 
raven colours indeed correlate with wsini. We investigate 
unevolved members of the Sco OB2 association in order 
to avoid the effects of stellar evolution on the colours. 
The oldest subgroup in Sco OB2 has an age of 14-15 Myr 
and this implies that the upper mass limit for stars still 
on the main sequence is about 12 M Q (using the models 
by Schaller et al.1992). This mass corresponds to stars 
of spectral type ~ Bl (see e.g., Harmanec 1988). In the 
following we concentrate on the B1-B9 stars to find out 
whether there are effects of rotation on Walraven colours. 

For the stars in our data we used the Walraven pho- 
tometry of de Geus et al. (1990) and the value of logT e ff 
and logg was taken from de Geus et al.(1989). For some 
stars we had to derive log T e g and log g from the photome- 
try and that was done as described in de Geus et al.(1989) 
and Brown et al. (1994). In Fig. 14a we show the HR- 
diagram in the logg-logT ff plane for all B1-B9 main se- 
quence dwarfs that are members of Sco OB2. Binaries were 
excluded from this diagram as well as B9p stars. The latter 
are known to have peculiar spectra and rotate slowly as a 
class. In Fig. 14a the symbol size scales with vsmi/v^. It 
is clear that for IogT c ff < 4.2 the stars with large values 
of fsim/fbr l ie above the main sequence. In this temper- 
ature range one finds B7-B9 stars. The Kelvin-Hclmholtz 
time-scale for these stars is ~ 2-5 x 10 5 yr. Even for mem- 
bers of US this is short enough for all B7-B9 stars to have 
evolved onto the zero-age main-sequence. Thus we do not 
expect any pre-main-sequence stars among this group. 

The only other reason for these stars to lie above the 
main sequence is that the transformation of the photome- 
try to logTeff and logg is in error. We therefore show the 
same sample of stars in the Walraven [B — U] vs. [B — L] 
diagram in Fig. 14b. The ZAMS is drawn in and the early 
type B-stars ([B — L] < 0.05) all lie on or very near this 
isochrone. In this colour-colour diagram the evolution of 



20 



A.G.A. Brown & W. Verschueren: Rotational Velocities of Early- Type Stars in Sco OB2 













0.1 































3 






(a) 













(b) - 




3.5 

M 

an 







I o° 


[B-U] 


0.1 
0.2 


- 
r 










A 






^ erg " 




0.3 






o\ 






4.5 






- 




0.4 


"I ■ ■ ■ ■ 1 . 












4.6 


4. 


1 4.2 4 




-0.05 


0.05 
LB 


0.1 
-L] 


0.15 0.2 



Fig. 14. (a) The Hertzsprung-Russell 
diagram in logT c ff and log(? for B1-B9 
dwarfs that are members of Sco OB2 
(b) The Walraven [B - U] vs. [B - L] 
colour-colour diagram. In both figures 
only single stars are included and the 
B9p stars are excluded. The symbol size 
scales with «sini/u br and isochrones 
(from Schaller et al.1992) are shown for 
0, 5, 10 and 15 Myr in (a) and for Myr 
in (b). 



a star away from the ZAMS is towards the lower left (see 
de Geus et al.1989). It is clear that for [B — L] > 0.05 
the stars are displaced away from the ZAMS in the direc- 
tion where evolved stars are expected. If we calculate the 
distance of each B7-B9 dwarf from the ZAMS (defined as 
the distance to the nearest point on the isochrone) and 
correlate that distance with vsini/vb r , we find a correla- 
tion at the 95.2% confidence level. Thus the photometry 
is indeed affected by rotation for the B7-B9 stars. 

The largest decrease in \ogg in Fig. 14a is 0.44 for a 
star with vsini/vb r = 0.84. The latter value implies that 
w > 0.96. Judging from Fig. 2 in Collins (1963) the value 
of the surface gravity then varies over a factor of ~ 3 from 
pole to equator. This is consistent with our observed factor 
of 2.8. We conclude that the Walraven colours are affected 
by stellar rotation. 

Note that some stars in Fig. 14b are located just above 
the ZAMS. It concerns a B2V, two B4V, a B8V and two 
B9V stars. One of the B4V stars is a star with emission 
lines. There is nothing peculiar about the other stars. The 
positions of the normal B2V and B4V stars can be ex- 
plained as observational errors. The B8 and B9 stars how- 
ever, lie further away from the ZAMS than the errors al- 
low. This can be explained by the fact that for B7-B9 
stars the empirical ZAMS for the Walraven system, as 
given in Brand & Wouterloot (1988), lies just above the 
ZAMS defined by Schaller et al.(1992). 

The lower surface gravity derived for the B7-B9 stars 
affects age determinations based on isochrone fitting as 
well as the derivation of mass distributions. The age of 
subgroups in the association will appear higher because of 
the apparent presence of evolved B7-B9 stars. However, 
we note that the ages derived by de Geus et al.(1989) are 
not affected by rotation. The stars in the B7-B9 range had 
little weight in the isochrone fitting procedure. The same 
holds for the ages derived for the Ori OBI subgroups by 
Brown et al.(1994), who excluded the B7-B9 stars from 
their isochrone fitting procedure. The IMFs derived by de 
Geus (1992) and Brown et al.(1994) for Sco OB2 and Ori 
OBI, respectively, are not affected by rotation. In both 



cases the mass range used by the authors excluded B7 B9 
stars. 



7. Conclusions and Future Work 

We have presented vsini values for 156 established or 
probable members of the Sco OB2 association. The pro- 
jected rotational velocities were derived from echelle spec- 
tra using three different techniques, depending on the an- 
ticipated value of vsini. We have shown that the results 
for the different techniques are consistent and that our re- 
sults are consistent with the standard system of SCBWP. 
On average the errors on our vsini values are ~ 10%. 

We analyzed our own data together with data for mem- 
bers of Sco OB2 taken from the literature. The analysis 
of the vsini distributions for the members of Sco OB2 
reveals that there are no significant differences between 
the subgroups in Sco OB2. The binary stars in Sco OB2 
rotate slower on the whole than the single stars. This is 
consistent with a picture in which angular momentum is 
transferred to the binary orbit during star formation, and 
possibly further removed due to tidal braking. For US we 
showed that the B7-B9 single stars rotate significantly 
faster than the B0-B6 stars. This may be consistent with 
the prediction about the distribution of rotational veloci- 
ties given by Mouschovias (1983). However, this result is 
derived from inhomogeneous data. 

Given the present data we cannot definitely exclude a 
form of random distribution of intrinsic rotational veloci- 
ties. The number of single stars in the data is not enough 
to get statistically significant results. Inclusion of the bina- 
ries complicates the picture due to the different physical 
mechanism by which binary components may acquire a 
low rotational velocity. 

The rectification of the vsini distributions by cither 
the technique of Lucy (1974, 1994) or that of Bcrnacca 
(1970), suggests that the peak in the distribution of rota- 
tional velocities shifts to higher values going from UCL to 
LCC to US. The interpretation of this finding in terms of 
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binary fraction per subgroup or evolutionary effects is not 
possible with the present data. 

For the Walraven photometric system we have now for 
the first time shown conclusively that the colours of stars 
are affected by rotation. This is most notably the case for 
the B7-B9 stars, where the distance to the ZAMS in the 
[£> — U] vs. [B — L] diagram correlates with vsini. However, 
we also conclude that the determination of ages and mass 
distributions for Sco OB2 (dc Geus et al.1989, de Geus 
1992) and Ori OBI (Brown et al.1994) are not affected by 
rotation. 

The first step in improving the analysis of the vsini 
distributions in Sco OB2 is to get a complete sample in 
all subgroups. This entails selecting members on the basis 
of kinematic data and performing a complete census of the 
binary population. It is very timely that the data from the 
HIPPARCOS astrometric mission will become available in 
early 1996. The SPECTER consortium in Leiden (see de 
Zeeuw et al.1994) was granted observing time on this mis- 
sion and we expect data for over 10 000 stars of spectral 
type earlier than F8 in the direction of nearby OB asso- 
ciations of which around 7000 are located in the direction 
of Sco OB2. The data will be analyzed for membership. 
The membership lists will greatly facilitate future exten- 
sion of the vsini data as well as searches for binaries. In 
addition to proper motions HIPPARCOS will also provide 
information on duplicity for binaries with separations of 
more than 0.1 arcsecond. 

On the theoretical side it will be very useful to calcu- 
late the expected effects of stellar rotation on Walraven 
colours. As a first step the calculations of CTC can be 
used. We have used the Roche model for rotating stars 
throughout this paper. This model is a very simple ap- 
proximation to reality and introduces a specific break- 
up velocity. In the case of non-uniformly rotating stars 
it is not clear that a break-up velocity exists and the ob- 
served equatorial velocity depends very much on the inter- 
nal distribution of angular momentum. For a given mass 
a stronger degree of concentration of angular momentum 
towards the centre of the star leads to a lower observed 
v c (see e.g., Bodenheimer 1971). So a consequence may 
be that there are many more low vsini stars than would 
be expected on the basis of uniformly rotating models. 
Furthermore, the displacement of stellar models in the 
HR-diagram changes and becomes more parallel to the 
main-sequence if the central concentration of angular mo- 
mentum increases. In highly differentially rotating models 
the stellar surface can develop a cusp at the poles and 
large variations in surface gravity may result. The present 
measurements of vsini still depend on the Roche model for 
rotating stars (for the rapid rotators) and on the correct- 
ness of approximating the atmosphere of a rotating star 
as plane-parallel pieces of Kurucz atmospheres on the sur- 
face of a rotating star. In order not to bias measurements 
of vsini better theoretical models are needed in which the 



atmospheric structure follows from the parameters of the 
rotating star. 

Finally, real constraints on the star formation process 
can only be obtained if these data can be compared to 
models that make detailed predictions of the distribution 
of rotational velocities, taking the star formation process 
into account. The models should take into account the 
early phases of star formation. In these phases the details 
of the (possible) magnetic braking process may determine 
whether a binary star is formed or not (e.g., Mouschovias 
1991). Subsequently the angular momentum history of the 
protostellar phase should be taken into account. Besides 
the angular momentum history, the origin of stellar masses 
and their radii should also be looked into, as these deter- 
mine the observed value of vsini. 
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Fig. 15. Frame of reference, attached to a rotating star, that 
is to be used in deriving the rotational broadening function. 
The Z-axis points towards the observer. 
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Appendix A: The Rotational Broadening Function 71 The R ota tional Broadening Function in In X-Space 



Following Tassoul (1978), we derive here the rotational 
broadening function A(x) for different limb-darkening 
laws. A rectangular frame of reference is used with its ori- 
gin at the centre of the star. The Z-axis points towards the 
observer, and the y-axis lies in the plane passing through 
the axis of rotation and the line of sight. It is assumed that 
the star is rotating uniformly, and that its observed disc 
is a circle of radius R (see also Fig. 15). For convenience: 
x = X/R and y = Y/R. 

According to assumption C of the classical model for 
rotating stars the line profiles in the spectrum of a star 
are constant over its surface and are given by: 



I(x, y;X- A c ) = r(A - X c )I Q (x, y) . 



(Al) 



where Iq (x, y) is the continuum intensity emergent at the 
stellar surface at position (x,y) and r(A — A c ) is the ra- 
tio of the intensity in the spectral line to the continuum 
intensity. The rotationally broadened spectral line is then 
defined as: 



+i 



v smz N 



f(A-A c )=/ r(A-A c -A c a: )A(x)dx. (A2) 

J-i c 

Here xvs'mi = —v z and A(x) is the rotational broadening 
function given by: 



A(x) 



Jo 



(l-x a ) 



2)1/2 



Io{x,y) dy 
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(l_ x 2)l/2 



I (x,y)dy 



(A3) 



To calculate A(x) a limb-darkening law has to be adopted. 
In Table 8 we list three forms for the limb-darkening law 
and the corresponding rotational broadening functions. In 
the first column the limb-darkening law is listed and in the 
second column the rotational broadening function. In the 
limb-darkening laws [i = cos 8 = yjl — (x 2 + y 2 ), where 
8 is the angular distance from the centre of the disc. The 
limb-darkening parameters are u, p, q, r, and s. 

The first two limb-darkening laws are the well known 
linear and quadratic approximation to limb-darkening and 
the third law is the square-root approximation proposed 
by Diaz-Cordoves & Gimcnez (1992). 

If 5(A) is the intrinsic spectrum, then the rotationally 
broadened spectrum S T (X) is given by: 



Rotational broadening of the spectrum depends on wave- 
length. This can be seen in the convolution integral in 
Eq. (A5) which is used to calculate 5 r (A). The wavelength 
dependence is in the sense that the rotational broadening 
function widens as A increases. This wavelength depen- 
dence is absent in In A-space (which is almost equivalent 
to pixel space for our spectra). By substitution of A' = In A 
into the convolution integral for 5 r (A) one finds: 



/■-ln(l-^fi) 

S;(ln A) = / 5' (In A — x')A'{x') dx' . 

J-ln(l+^iiH) 



(A6) 



where S"(lnA) is the spectrum and A'(x') the rotational 
broadening function in In A-space. 

The rotational broadening functions in In A-space for 
the different limb darkening laws can be found from the 
rotational broadening functions in A-space (listed above) 
as follows: 



.1':,'. —A 1 ' 



vsini/c \ vsini/c 
vsini 



-In II 



, , / _ vsmi , 
) < x < <-^[l- — \ ,A,) 

This can be derived from the following: 



S r (e V ) = ^(A') = 
Setting u = z/e x , one finds: 



S(e x '-z)^-A ( f - I <fo.(A8) 
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(A9) 



since S(y) — S'(lny). Substituting x' = — ln(l — u) finally 
yields: 



-ln(l-^a-) 




vsini/c \ vsini/c 



dx' . 



(A10) 



Sr(X) 



S X 



, vsmi . , 
Xx I A(x) dx , 



or, substituting z — Xx^^-: 



Sr(X) -- 

where A r = A^. 



S (A-z)lA(f ) d,. 



( A4) ^- ^- The Width of Rotationally Broadened Lines 

At low values of the projected rotational velocity (t;sim < 
80 kms -1 ) the measured width of a line that is not af- 
fected by other broadening mechanisms (Apart from ther- 
. . mal and instrumental) can be directly converted into a 
^ ' value of wsini. In this paper we use metal lines for this 
purpose and measure their FWHM by fitting Gaussians. 
This implies a measurement of the second moment of the 
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line profile. If the second moment of a Gaussian is a then 
the FWHM is given by 2\/21n2cr. The second moment of 
a broadened, intrinsically infinitely sharp line centered on 
Ao is given by: 

rA +A t 



<r 2 (A ) = / ° ' (A — Ao) 2 f(A) dX— 

JXo-\t 



r A +A r 



(A - Ao)f(A) dX , (All) 
with f (A) in this case given by: 

f(A) = j + * 5(A - A - dz . (A12) 



Because of symmetry Eq. (All) reduces to: 
f +i 

cr 2 (A ) = J X 2 x 2 A(x)dx. 



(A13) 



Note that the last integral contains an integrand in which 
only terms of the form x 2 {l — x 2 ) u occur, which are all 
symmetrical around zero. The Beta-function can be given 
as: 

B(x,y) = 2 [ t^-^l-ty^dt 
Jo 

(Rc(x) > 0, Rc(y) > 0) . (A14) 

So one can easily evaluate a 2 using this expression for the 
Beta-function and the relation between the Beta function 
and the Gamma-function. For linear limb-darkening the 
result is: 



( (A r = Ao=). (A15) 



The results for quadratic and square-root limb-darkening 
are: 



a 2 (A ) = A r 2 



a 2 (A ) = A r 2 



10 



P 



6-2p-q 

5 7^ 13 



15 - 5r - 3s 



(A16) 
(A17) 



The ECHELEC spectra that are being analyzed are 
recorded in In A-space and in that case one can show that 
the general equation for the second moment a' is: 



l\2 



+1 , 2/-, wsini 

ln 2 (l y)A(y)dy- 

-l c 

+i • • \ 2 

ln(l y)A{y)dy 

-l c 



(A18) 



The two integrals in the equation cannot be calculated 
straightforwardly. However, in the limit vsim/c -C 1 the 
integrals reduce to the same form as for a but without 
the dependence on wavelength. Since the method is specif- 
ically applied in cases of low values of vsmi the approxima- 
tion should hold when determining t>sim from the FWHM 
of the metal lines. For B3-B5 stars (T off w 15 000K) the 
values of u, p, q, r, and s in the Stromgren w-band at 
\ogg = 4.5 are 0.408, 0.159, 0.324, -0.102, and 0.723, 
respectively (taken from Diaz-Cordoves et al.1995). This 
leads to the following values for the coefficients of A r in 
Eqs. (A15), (A16) and (A17); 0.234, 0.210, and 0.236. This 
illustrates that the derived values of vsmi are not very 
sensitive to the choice of limb-darkening law. 

Appendix B: Fourier Filtering of Template and 
Target Spectrum 

As mentioned in Sect. 3.2 it is necessary to filter both 
the target and the broadened template spectrum before 
comparing the two. The filtering is done in Fourier-space 
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and the high-pass filter is constructed as follows. The in- 
formation relevant to the rotation of a star can be found 
in its spectrum on spatial scales ~ 2vsini (in kms -1 ). 
All components in the spectrum with longer spatial scales 
than this value carry essentially no information on vsini 
and can be filtered out. So based on the highest value in 
the vsini search interval one can construct the high-pass 
filter in Fourier-space. This filter is a step function with 
its edge tapered by a narrow Gaussian. The edge of the 
filter is placed at the frequency corresponding to the up- 
per limit of the tested vsini values. In practice two extra 
frequencies below the edge value are retained. This is to 
insure that one does not throw away too much informa- 
tion. After filtering, the comparison between template and 
target is done as described in Sect. 3.2. 

The data do not contain a number of points equal to 
a power of two as required by the FFT routine used. So 
in practice in order to avoid creating large discontinuities 
in the data by zero padding we connect the endpoints in 
the data by a straight line and subtract that line from the 
data. The endpoints are often very noisy so we first dis- 
card the outer 10 points at each end. Components with a 
spatial wavelength smaller than ~ 260 kms -1 were never 
filtered out in order to avoid being dominated by small 
mismatches between template and target spectrum. Sig- 
nal to noise considerations also lead us to throw out 40 
extra pixels at both edges of every spectral order (con- 
taining 500 pixels in total) before calculating the mean 
absolute difference between template and target. 

Note that before doing the comparison one has to make 
sure that the template and target spectrum are shifted to 
the same radial velocity. It turns out that aligning the 
spectra to within one pixel is sufficiently accurate (see 
Appendix C). This can be done by eye. 

Appendix C: The Effects of Gravity- and Limb- 
Darkening on Template Broadening 

As discussed in Sect. 3.2 we expect gravity- and limb- 
darkening to have more severe effects when one uses tem- 
plate broadening to derive vsini. We investigated this by 
generating synthetic data from the ab initio models by 
CTC. We took the models of the Hei A4026 line and the 
H7 line and constructed synthetic 'spectral orders' with 
the same characteristics as an ECHELEC spectral order. 
The continuum was taken to be constant (normalized flux 
1) and Poisson noise was added artificially. The artificial 
spectral order is calculated in In A-space because the mod- 
els are given as a function of wavelength. Models with low 
values of vsini were used as templates. This investigation 
was done for B3 and B9 models in order to compare results 
for different spectral types. 

An example is shown in Fig. 16 for the He I A4026 line. 
In both cases the target spectrum is that of a star with 
w = 0.8 and i = 90° The corresponding values of vsini are 



250 kms _1 and 213 kms -1 , respectively, for the B3 and B9 
target. The noise in the target spectra corresponds to a 
signal to noise ratio of 250. This is a typical value for or- 
der 142. The templates were non-rotating models in both 
cases. The template broadening technique applied in these 
two cases yields values for vsini of 240 kms _1 for the B3 
target and 180 kms _1 for the B9 target. The broadened 
template spectra in Fig. 16 correspond to these values. 
Also plotted in Fig. 16 are the same model spectra without 
the noise in order to show the actual differences between 
the broadened template and the target. In both cases one 
can see that the He 1 A4026 line from the ab initio model is 
narrower and less deep than the corresponding classically 
broadened template. This effect is largest for the B9 star. 
In both cases vsini is underestimated by template broad- 
ening, but the error for the B3 star is less than 5%. Note 
that for both spectral types the fit of the noisy broadened 
template to the noisy data looks good in the high-pass 
filter case even though the value of vsini is wrong. 

For the parameters of the example shown (S/N = 250, 
no radial velocity shifts, zero- rotation template), template 
broadening gives the correct results up to w = 0.8 for 
B3 stars (vsini is underestimated by ~ 5% for w = 0.8, 
i = 90°). For B3 stars rotating at w > 0.9 the value 
of vsini can be underestimated up to 15%. For B9 stars 
template broadening only works well up to w = 0.5 (vsini 
is underestimated by ~ 10% for 10 = 0.5, i = 90°). For 
higher values of w, vsini is underestimated by at least 10- 
15% and for w > 0.9 the discrepancies between template 
broadening and model fitting are so severe that depending 
on the inclination vsini can be overestimated as well as 
underestimated by ~ 20%. 

For the B3 models we investigated the effects of other 
parameters in the synthetic data. At a lower signal-to- 
noise ratio in the spectrum (S/N=100) the results are not 
affected. If one introduces small radial velocity shifts of 1 
to 2 pixels between template and object the vsini values 
from template broadening are off by 5-8% at the lowest 
intrinsic vsini values. At higher values of vsini the results 
are not affected; it becomes less important at those values 
to have the phases of template and object match exactly. 
If one uses a template for which vsini w 30 kms -1 the 
value for the target is underestimated even more. For w < 
0.8 vsini is now underestimated by 5-10% and for higher 
values of w the underestimates can be up to 17%. 

The results can be summarized by stating that the 
vsini values for template broadening of the Hei A4026 
line can be trusted up to ~ 200 kms -1 for early type B3 
stars. For the B9 stars one should only trust template 
broadening for values of vsini less than 120-150 kms -1 . 

The situation is reversed when using the H7 line. Now 
it is the B3 stars for which vsini is underestimated at low 
values of w. The value of vsini from template broadening 
can be off by as much as 30% at w = 0.5. For B9 stars the 
value of vsini from template broadening can be trusted 
up to high values of w. 
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Fig. 16. Example of template broadening applied to synthetic data. The four panels on the left indicate the results for B3 
models. The panels on the right show the results for B9 models. The top four panels show a direct comparison between the 
best fitting broadened template (thick line) and target (thin line) . The best fits are based on the filtered data. In the top two 
panels the dashed line is the broadened template corresponding to the real value of -usim of the object. The bottom panels show 
the comparison after high-pass filtering of the data. The same results are shown for synthetic data with and without noise. 
The noisy data have a signal to noise ratio of 250. In both cases the template is a non-rotating star and the target rotates at 
80% of the critical angular velocity. The inclination angle is 90° in both cases. The resulting values of vsim are 250 kms _1 and 
213 kms -1 , respectively, for the B3 and B9 targets. The values of usim found from template broadening are 240 kms" 1 and 
180 kms -1 , respectively. The broadened templates in this figure correspond to the latter vsini values. Note the differences in 
vertical scale. 



The difference in behaviour of He I A4026 and H7 with 
spectral type is due to a difference in sensitivity to surface 
gravity. The Balmcr lines are most sensitive to surface 
gravity changes for the early type B-stars. A reduction of 
the surface gravity causes the wings of the Balmcr lines 
to become less important whereas the line strength in the 
core remains the same. This effect leads to underestimated 
vsini values and is strongest for the early-type B-stars. 
The He 1 A4026 line shows the opposite behaviour and it 
is the late-type B-stars for which the line is most sensitive 



to changes in surface gravity. In this case the whole line 
becomes weaker and narrower. 

The general conclusion is that one should be careful 
when using template broadening results as soon as wsinz 
is larger than ~ 120 150 kms -1 . Visual inspection of the 
results is required and one should take the spectral type 
of the target star into account. 

The role of limb-darkening was already discussed in 
Sect. 3.1. In the simulations we also used the linear limb- 
darkening law with parameter u = 0.4 to calculate the 
broadened template. The results suggest that the choice 
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of limb-darkening law does not affect the outcome by 
much. To investigate this a little further we re-determined 
vsini for a couple of our stars with a different choice of 
limb-darkening law. The two other forms we tried are 
the quadratic and square-root approximations to limb- 
darkening (see e.g., Dfaz-Cordoves & Gimcnez 1992). We 
took a couple of stars and used a limb-darkening law and 
parameters appropriate for the spectral type of the stars. 
The results were not affected and we are confident that 
the choice of limb-darkening law is not important when 
using template broadening. The reason is that the differ- 
ences due to limb-darkening are obscured by the noise in 
the spectra. 

Appendix D: Notes to Table 2 

HD 94650: The Hei A4026 line of this star is best fit with 
a B3 model. The appearance of the Helium lines in order 
142 also points to an earlier spectral type. 

HD 95122: The B5 model for Hei A4026 gives the best fit. 

HD 100262: This star is a supergiant and the measured 
vsini may actually be determined by macro-turbulence in 
the atmosphere of this star. 

HD 101947: The value of vsini for this star was estimated 
by eye from its spectrum. This value may also be a mea- 
sure of macro-turbulence rather than rotation. Note that 
this star is also classified as G01a + (J. Lub, personal com- 
munication) . 

HD 104878: For this star B9 models were used for 
Hei A4026. 

HD 105382: The best template for the template broaden- 
ing method was HD 120709, which is a B5 star. Indeed, the 
star is classified as B6IIIe in the Bright Star Catalogue. 

HD 105580: The best fitting template is HD 132955, a 
B3V star. The spectrum of this star is clearly not of type 
B6. Houk and Cowley (1975) classify this star as B3IV. 

HD 107566: The value of vsini is taken from Uesugi & 
Fukuda (1981) and revised downwards by 5% to make 
the value consistent with the standard system given by 
SCBWP. 

HD 107696: The best fitting model for Hei A4026 is for 
spectral type B5. The strength of the Hei A4009 line is 
clearly not consistent with spectral type B9. Houk & Cow- 
ley (1975) list B7Vn as spectral type for this star. 

HD 108248: HD 149438 (B0V, r Sco) was used as tem- 
plate. 

HD 108249: This star is a double-lined spectroscopic bi- 
nary. The model fits were done by direct comparison 
(without Fourier filtering) of the strongest component to 
the models. This makes this vsini value rather uncertain. 

HD 108257: The Hei A4026 line in this star is broader than 
the models. The value of vsini may thus be overestimated. 



HD 108483: The best fitting template is HD 96706, a B2V 
star. The spectral type listed in Houk (1978) is B2V. 

HD 109787: The value of vsini is taken from SCBWP. 

HD 110432: The best fit of the Hei A4026 line is for a B5 
model. 

HD 111123: HD 149438 was used as template for method 
2. HD 11123 is a non-radial pulsator and a (3 Cephei vari- 
able. 

HD 112078: B3 models were used to fit the Hei A4026 
line. The line profile may be influenced by the presence of 
a secondary in the spectrum. 

HD 112091: B4 models were used to fit the Hei A4026 line. 

HD 112409: B7 models were used to fit Hei A4026. 

HD 113703: B3 models were used to fit Hei A4026. Houk 
(1978) lists spectral type B4V for this star. 

HD 113791: Listed as B3V by Houk (1978). 

HD 113902: B7 models were used to fit Hei A4026. 

HD 114529: The value of vsini for this star was estimated 
by eye. It may be a double-lined spectroscopic binary. 

HD 116072: B3 models were used to fit Hei A4026. The 
spectral type of this star is B2.5Vn according to the Bright 
Star Catalogue. 

HD 116226: B5 models were used to fit Hei A4026. Houk 
(1978) gives B5V for the spectral type of this star. 

HD 118716: B3III models were used for method 3. This 
star is a non-radial pulsator and a (3 Cephei variable. 

HD 118978: B8 models were used to fit Hei A4026. 

HD 119361: B7III models were used to fit Hei A4026. 

HD 120324: This star is a non-radial pulsator. 

HD 121190: B7 models were used to fit Hei A4026. 

HD 121263: The value of vsini may be influenced by the 
presence of a secondary component in the spectrum. 

HD 124367: B3 models were used to fit Hei A4026. 

HD 125238: The stellar Hei A4026 line is much deeper 
than the B3 model line. The value of vsini may thus be 
underestimated. 

HD 125745: B6 models were used to fit Hei A4026. The ap- 
pearance of the Helium lines in order 142 clearly suggests 
an earlier spectral type. 

HD 126341: This star is a non-radial pulsator and a (3 
Cephei variable. 

HD 126981: B6 models were used to fit Hei A4026. The 
appearance of the Helium lines in order 142 suggests an 
earlier spectral type. 

HD 127381: HD 126341 (B2IV) was used as template. 

HD 127972: The value of vsini is taken from SCBWP. 

HD 128345: B4 models were used to fit Hei A4026. Houk 
(1978) lists the spectral type of this star as B3/B4V 

HD 129056: This star is (3 Cephei variable. 
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HD 131625: B9 models were used to lit Hei A4026. The 
models do not fit well. 

HD 132058: HD 126341 was used as template. Houk (1978) 
lists B2IV as spectral type for HD 132058. 

HD 132851: The value of wsinz was estimated from a visual 
inspection of the spectrum. 

HD 133937: B3 models were used to fit Hei A4026. The ap- 
pearance of the Helium lines in order 142 clearly suggests 
an earlier spectral type. Houk (1978) lists the spectral type 
as B5/B7V. 

HD 135876: B6 models were used to fit Hei A4026. The ap- 
pearance of the Helium lines in order 142 clearly suggests 
an earlier spectral type. Houk (1978) lists the spectral type 
as B7V. 

HD 136298: B3 models were used to fit Hei A4026. This 
star is a f3 Cephei variable. 

HD 136664: B3 models were used to fit Hei A4026. 

HD 136933: The value of vsmi was estimated from a visual 
inspection of the spectrum. 

HD 137058: B9 models were used to fit Hei A4026. The 
models do not fit well. 

HD 137432: HD 132955 (B3V) was used as template. The 
Bright Star Catalogue lists the spectral type of this star 
as B4Vp. 

HD 140784: B4 models were used to fit Hei A4026. The 
appearance of the Helium lines in order 142 clearly sug- 
gests an earlier spectral type. The Bright Star Catalogue 
lists the spectral type as B7Vn. 

HD 141637: B3 models were used to fit Hei A4026. The 
spectrum shows signs of a secondary component. 

HD 142114: The stellar Hei A4026 line is much deeper 
than the B3 model line. The value of vs'mi may thus be 
underestimated. 

HD 142184: The stellar Hei A4026 line is much deeper 
than the B3 model line. The value of vsmi may thus be 
underestimated . 

HD 142629: The value of vs'mi is taken from Uesugi & 
Fukuda (1981) and revised downwards by 5% to make 
the value consistent with the standard system given by 
SCBWP. 

HD 142983: The value of vsmi was estimated from a visual 
inspection of the spectrum. According to SCBWP this star 
has a projected rotational velocity of 400 km s^ 1 . However, 
judging from the spectrum this clearly is not the case. This 
star is also a non-radial pulsator. 

HD 143018: The value of vsmi is taken from the Bright 
Star Catalogue and the spectrum was inspected visually 
to check whether the value listed is reasonable. 

HD 143275: Bl models were used to fit Hei A4026 and HD 
149438 was used as template for method 2. The Bl model 
profiles are not really deep enough to fit the data. 

HD 144217: HD 149438 was used as template. 



HD 144294: The stellar Hei A4026 line is much deeper 
than the B3 model line. The value of vsmi may thus be 
underestimated. 

HD 144987: B7 models were used to fit Hei A4026. 

HD 145482: The value of vsmi is influenced by the fact 
that this star is a double-lined spectroscopic binary. 

HD 147084: The value of vsmi was estimated from a visual 
inspection of the spectrum. 

HD 147165: This star is a (3 Cephei variable. 

HD 147628: B6 models were used to fit Hei A4026. The 
appearance of the Helium lines in order 142 suggests an 
earlier spectral type. 

HD 147933: The value of vsmi may be influenced by the 
presence of a secondary component in the spectrum. 

HD 148184: B4 models were used to fit Hei A4026. This 
star shows emission cores in the Balmer lines present in 
our spectra. 

HD 149757: The value of vsmi is taken from SCBWP. 
This star is a non-radial pulsator. 

HD 151804: The value of vsmi was estimated from a visual 
inspection of the spectrum. 

HD 151890: Bl models were used to fit Hei A4026. The 
value of vsmi may be influenced by the secondary com- 
ponent in the spectrum. 

HD 154204: B4 models were used to fit Hei A4026. The 
appearance of the Helium lines in order 142 suggests an 
earlier spectral type. The Bright Star Catalogue lists the 
spectral type as B6IV. 

HD 157042: The value of vsmi was estimated from a vi- 
sual inspection of the spectrum. This star is a non-radial 
pulsator. 

HD 157056: This star is a (3 Cephei variable. 

HD 157246: The value of vsmi is taken from SCBWP. 
This star is a non-radial pulsator. 

HD 158427: The value of vsmi may be influenced by the 
presence of a secondary component in the spectrum. 

HD 158926: The value of vsmi was derived from the 
spectrum by visual inspection. The vsim value listed in 
SCBWP is 190 kms- 1 . This star is a Cephei variable. 

HD 168905: The stellar Hei A4026 line is much deeper 
than the B3 model line. The value of vsmi may thus be 
underestimated. 

HD 171034: HD 126341 was used as template. 
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Table 2a. Projected rotational velocities in Lower Centaurus Crux. 



HD 


HIC 


MK 


v sin i 


<t(v sin?) 


^br 


v sin % 


Meth. 


Mem. 


Dupl. 










kms —1 


kms —1 


kms -1 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




94650 


53272 


B6V 


400 


13 


411 


0.973 


3 








95122 


53556 


B7V 


152 


16 


388 


0.392 


3 




RV 




96706 


54327 


B2V 


44 


2 


459 


024 


4 




CON 




99103 


55597 


B5 


102 


4 


388 


0.263 


3 








99264 


55657 


B2.5IV-V 


125 


5 


411 


0.304 


3 




RV 




99556 


55831 


B3IV 


188 


10 


411 


0.457 


3 




CON 




100262 


56250 


A3Ia 


30 


3 


152 


0.197 


1 










56561 


B9II* 


155 


17 


278 


558 


3 








100929 


56606 


B2.5IV 


23 


4 


459 


0.050 


1 




CON 




101947 


57175 


F9Ia 


30 


10 


122 


0.246 


4 






* 


102776 


57669 


B3V 


251 


9 


411 


0.611 


3 


m 


RV 




103079 


57851 


B4V 


47 


7 


401 


0.117 


1,2 


M 






iUOOO^l 


OooZD 


DO V 


1 A Q 


o 
o 


/ill 


U.ODO 


9 ^ 


1V1 


POM 




103961 


58379 


B8III 


60 


7 


288 


0.208 


1,2 








104841 


58867 


B2IV 


25 


2 


459 


0.054 


1 


m 


SB1 




104878 


58884 


A0V 


225 


30 


350 


0.643 


3 




RV 




105382 


59173 


B2IIIne 


75 


13 


303 


0.248 


1,2 


M 


RV 






59196 


DZ1 V lie 


298 


14 


411 


725 


3 


M 


RV 




105580 


59266 


B6V 


44 


8 


411 


0.108 


1,2 






* 


105937 


59449 


B3V 


129 


6 


411 


0.314 


3 


M 


RV 




106490 


59747 


B2IV 


135 


6 


411 


0.328 


3 


M 


CON 




106983 


60009 


B2.5V 


65 


9 


459 


0.142 


1,2 


M 






1U 1 ODD 


DUOZU 


Am 


114 


44 


357 


319 








* 


107696 


60379 


B9V 


382 


40 


388 


0.985 


3 


m 




* 


108248 


60718 


B0.5IV 


88 


5 


504 


0.175 


2 




SB1 


* 


108249 




B1V 


150 


20 


491 


0.305 


3 




RV 




108250 




B4V 


131 


6 


401 


0.327 


3 




SB1 




1 UoZ O i 


60710 


B3Vn 


298 


12 


411 


725 


3 


M 


CON 


* 


108483 


60823 


B3V 


169 


12 


459 


0.368 


2 


M 


RV 


* 


109026 


61199 


B5V 


188 


10 


401 


0.469 


3 


M 


RV 




109668 


61585 


B2IV-V 


114 


11 


459 


0.248 


2 


M 


RV 




109787 


61622 


A2V 


280 


28 


365 


0.767 


4 








110304 


61932 


A1IV 


79 


3 


337 


0.234 


2 




RV 




110432 


62027 


B2pe 


388 


42 


388 


1.000 


3 




RV 


* 


110879 


62322 


B2V 


139 


6 


411 


0.338 


3 


m 


RV 




110956 


62327 


B3V 


26 


4 


411 


0.063 


1 


M 


RV 




111123 


62434 


B0.5III 


40 


5 


401 


0.100 


12 


M 


RV 


* 


112078 


63007 


B4Vn 


298 


27 


411 


0.725 


3 


M 


RV 


* 


112091 


63005 


B5Vne 


242 


10 


401 


0.603 


3 


M 




* 


112092 


63003 


B2IV-V 


34 


3 


459 


0.074 


1 


M 


CON 





Table 2a. (Continued) 



HD 


HIC 


MK 


v sin i 


<t(v sin?) 


^br 


v sin % 


Meth. 


Mem. 


Dupl. 










kms -1 


kms -1 


kms -1 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




112409 


63210 


B8V 


164 


25 


373 


0.440 


3 








113703 


63945 


B5V 


140 


7 


411 


0.341 


3 


M 






113791 


64004 


JdI.o V 


15 


2 


475 


0.032 


1 


1V1 


C ID 1 


* 


113902 


64053 


B8V 


281 


46 


373 


0.753 


3 






* 


114529 


64425 


B8V 


250 


25 


363 


0.689 


4 




RV 


* 


115823 


65112 


B6V 


42 


5 


381 


0.110 


1 


M 


RV 




115846 


65181 


B3V 


168 


7 


411 


0.409 


3 








116072 




B4V(V?) 


223 


11 


411 


0.543 


3 






* 


116087 


65271 


B3V 


223 


11 


411 


0.543 


3 


M 


CON 




116226 


65303 


B6IV 


215 


22 


388 


0.554 


3 






* 


118716 


66657 


B1III 


114 


15 


377 


0.302 


2,3 


M 


RV 


* 


118978 


66849 


B9IV 


182 


29 


363 


0.501 


3 


m 




* 


124471 


69763 


B1.5III 


38 


4 


377 


0.101 


1 




RV 





Table 2b. Projected rotational velocities in Sco OB2-1. 



HD 


HIC 


MK 


v sin i 


<t(v sin?) 


^br 


v sin % 
^br 


Meth. 


Mem. 


Dupl. 










kms -1 


kms -1 


kms -1 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




151804 


82493 


09e 


75 


15 


508 


0.148 


4 






* 


152234 


82676 


B0.5Ia 


90 


11 


400 


0.225 


1 




RV 




152236 


82671 


Bllae 


75 


2 


373 


0.201 


1 








157042 


85079 


B2IIIne 


275 


30 


364 


0.755 


4 






* 


157246 


85267 


Bllb 


230 


23 


382 


0.602 


4 




RV 


* 


158408 


85696 


B2IV 


28 


4 


459 


0.061 


1 




RV 




158427 


85792 


B2Vne 


295 


19 


411 


0.718 


3 




RV 


* 


158926 


85927 


B1.5IV+... 


150 


20 


475 


0.316 


4 




SB2 


* 


163685 


88012 


B3II/III 


13 


3 


326 


0.040 


1 








165024 


88714 


B2Ib 


108 


3 


329 


0.328 


1 




RV 




166596 


89290 


B2.5III 


212 


9 


326 


0.650 


3 








167128 


89605 


B3IIIpe 


47 


6 


326 


0.144 


1 




RV 




167756 


89726 


B0.5Ia 


76 


6 


400 


0.190 


1 








168905 


90200 


B2.5Vn 


248 


12 


411 


0.603 


3 




CON 


* 


169467 


90422 


B3IV 


23 


3 


411 


0.056 


1 




RV 




170523 


90853 


B3III 


11 


2 


326 


0.034 


1 




SB1 




171034 


91014 


B2III/IV 


107 


10 


459 


0.233 


2 




CON 


* 


172910 


91918 


B2V 


12 


2 


459 


0.026 


1 




RV 





Table 2c. Projected rotational velocities in Upper Centaurus Lupus. 



HD 


HIC 


MK 


v sin i 


<t(v sin?) 


^br 


v sin % 


Meth. 


Mem. 


Dupl. 










kms —1 


kms —1 


kms -1 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




119361 


66984 


B8III 


139 


18 


296 


0.470 


3 








120307 


67464 


B2IV 


65 


6 


459 


0.142 


1,2 


M 


SB1 






67472 


R9TV 

13 Zl V - Vc 


130 


g 


459 


283 


2 3 


M 


RV 


* 


120640 


67663 


B2Vp 


21 


3 


459 


0.046 


1 




CON 




120709 


67669 


B5 


8 


3 


388 


0.021 


1 


M 


RV 




120710 




B9V 


175 


28 


350 


0.500 


3 


m 


SB1 




120955 


67786 


B4IV 


40 


2 


401 


0.100 


1 


M 


SB1 




121190 


67973 


B8V 


240 


32 


373 


643 


3 






* 


121263 


68002 


B2.5IV 


128 


6 


411 


0.311 


3 




SB2 


* 


121743 


68245 


B2IV 


79 


7 


459 


0.172 


2 


M 


RV 




121790 


68282 


B2IV-V 


124 


7 


459 


0.270 


2 


M 


CON 




122980 


68862 


B2V 


15 


3 


459 


0.033 


1 


M 


RV 




1 Z^loO / 


oyoio 


d*± v lie 


^9^ 
OZO 


ia 
o^l 


/ill 
^ii i 


U. / OD 


O 


in 


POM 


* 


125238 


69996 


B2.5V 


222 


11 


411 


0.540 


3 


in 


CON 


* 


125745 


70243 


B8V 


271 


26 


381 


0.711 


3 






* 


125823 


70300 


B2V 


15 


2 


459 


0.033 


1 


M 


RV 




126341 


70574 


B2IV 


15 


2 


459 


0.033 


1 




RV 






70915 


B8Vn 


320 


40 


381 


840 


3 






* 


127381 


71121 


B2III 


69 


10 


364 


0.190 


2 






* 


127971 


71353 


B6V 


177 


21 


381 


0.465 


3 








127972 


71352 


BlVn+... 


260 


26 


491 


0.530 


4 


M 


RV 




128345 


71536 


B5V 


186 


11 


401 


0.464 


3 




CON 




1 900^ 


71860 


Bl 5III 


16 


3 


377 


042 


^ 


M 


RV 


* 


129116 


71865 


B2.5V 


129 


6 


411 


0.314 


3 


M 


CON 




130807 


72683 


B5IV 


27 


3 


388 


0.070 


1 


M 


RV 




131120 


72800 


Dktt /TTT 

B711/111 


57 


6 


296 


0.193 


1 


m 


CON 




131562 


73095 


A2III 


57 


3 


307 


0.186 


2 




RV 




1 ^1 ft9 ^ 

1 1 OZ O 


73049 


A0V 


183 


12 


350 


523 


3 




SB2 


* 


132058 


73273 


B2III 


92 


6 


364 


0.253 


2 


M 


RV 


* 


132200 


73334 


B2IV 


32 


3 


459 


0.070 


1 


M 


RV 




132851 


73566 


A4IV 


120 


18 


350 


0.343 


4 








132955 


73624 


B3V 


8 


1 


411 


0.019 


1 


M 


RV 




1 ^^Q 1 ^? 

looyo / 


7 a i nn 

/ ^IIUU 


R7V 
n / V 


A (19 
^1UZ 


i ^ 
lo 


/ill 
^ii i 


n Q7S 

u.y / o 


O 


1V1 


RV 


* 


133955 


74117 


B3V 


135 


14 


411 


0.328 


3 


M 


RV 




134481 


74376 


B9V 


175 


28 


350 


0.500 


3 








134687 


74449 


B3IV 


13 


3 


411 


0.032 


1 


m 


SBl 




looo / C) 


7/1 q t^n 


D3 V 


lie 
llo 


y 


QC1 


n qi n 


Q 
O 




CDO 
ODZ 


* 


136298 


75141 


B1.5IV 


193 


9 


411 


0.470 


3 


M 


RV 


* 


136504 


75264 


B2IV-V 


41 


6 


459 


0.089 


1 


m 


SB2 




136664 


75304 


B4V 


177 


17 


411 


0.431 


3 


M 


RV 


* 


136933 


75439 


Ap... 


35 


5 


362 


0.097 


4 






* 


137058 


75501 


A0V 


268 


58 


350 


0.766 


3 






* 


137432 


75647 


B5V 


77 


12 


411 


0.187 


2 


M 


SBl 


* 


138769 


76371 


B3IVp 


67 


6 


411 


0.163 


1 


M 


RV 





Table 2d. Projected rotational velocities in Upper Scorpius. 



HD 


HIC 


MK 


v sin i 


<t(v sin?) 


^br 


v sin % 
fbr 


Meth. 


Mem. 


Dupl. 










kms -1 


kms -1 


kms -1 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




139365 


76600 


B2.5V 


134 


8 


459 


0.292 


2 


M 


SB2 




140008 


76945 


B5V 


11 


3 


388 


0.028 


1 


M 


SB2 






77286 


B8V 


401 


40 


401 


1 000 


3 






* 


141556 


77634 


B9.5III-IV 


8 


2 


325 


0.025 


1 




SB2 




141637 


77635 


B1.5Vn 


227 


12 


411 


0.552 


3 


M 


RV 


* 


142096 


77811 


B3V 


146 


5 


411 


0.355 


3 


M 


SB2 




142114 


77840 


B2.5Vn 


240 


12 


411 


0.584 


3 


M 


CON 




142184 


77859 


B2V 


255 


17 


411 


620 


3 


M 


CON 


* 


142629 


78105 


A3V 


81 


10 


368 


0.220 


4 






* 


142669 


78104 


B2IV/V 


98 


8 


459 


0.214 


2 


M 


SBl 




142883 


78168 


B3V 


14 


2 


411 


0.046 


1 


M 


SB2 




142983 


78207 


B8Ia/Iab 


40 


10 


187 


0.214 


4 


M 






1 A 9QQH 

i^izyyu 
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DO V 


1 78 
1 / o 


99 
ZZ 


OOO 


n a t^Q 


O 


1V1 


RV 




143018 


78265 


B1V+... 


100 


15 


491 


0.204 


4 


M 


SB2 


* 


143118 


78384 


B2.5IV 


191 


8 


411 


0.465 


3 


M 


RV 




143275 


78401 


B0.2IV 


148 


8 


489 


0.303 


2,3 


M 


SBl 




144217 


78820 


B0.5IV 


91 


8 


504 


0.181 


2 


M 


SB2 




144218 


78821 


B2V 


56 


g 


459 


122 


1 2 


M 


RV 




144294 


78918 


B2.5Vn 


252 


7 


411 


0.613 


3 


M 


RV 


* 


144470 


78933 


B1V 


100 


6 


491 


0.204 


2 


M 


CON 




144844 


79098 


B9V 


10 


2 


350 


0.029 


1 


M 


SB2 




144987 


79199 


B8V 


240 


32 


373 


0.643 


3 




RV 




145482 


79404 


B2V 


174 


10 


459 


379 


2 


M 


SB2 


* 


145483 


79399 


B9V 


225 


30 


350 


0.643 


3 


M 






145502 


79374 


B2IV 


162 


8 


411 


0.394 


3 


M 


SBl 




146624 


79881 


A0V: 


41 


3 


362 


0.113 


1,2 








147084 


80079 


A4II/III 


10 


5 


261 


0.038 


4 


M 


RV 




147165 


80112 


B1III 


56 


15 


390 


144 




M 


SB2 


* 


147628 


80390 


B8V 


170 


12 


381 


0.446 


3 




RV 


* 


147933 


80473 


B2V 


196 


10 


411 


0.477 


3 


M 


RV 


* 


147934 




B2V 


223 


11 


411 


0.543 


3 


M 


RV 




148184 


80569 


B2Vne 


148 


15 


401 


0.369 


3 


M 


SBl 




148703 


80911 


R9TTT TV 

13 Zlll-1 V 


70 


g 


364 


192 


2 


M 


RV 




149438 


81266 


B0V 


10 


2 


515 


0.019 


1 


M 


RV 




149757 


81377 


09.5V 


320 


32 


601 


0.532 


4 


M 


CON 


* 


151890 


82514 


B1.5IV+... 


180 


9 


411 


0.438 


3 


M 


SB2 


* 


151985 


82545 


B2IV 


52 


5 


459 


0.113 


1,2 


M 


RV 




154090 


83574 


B2Iab 


78 


2 


314 


0.248 


1 








154204 


83567 


B7IV/V 


372 


24 


401 


0.928 


3 






* 


155889 


84444 


Bl/B2Ib/II 


38 


4 


361 


0.105 


1 








157056 


84970 


B2IV 


31 


3 


459 


0.068 


1 


M 


RV 


* 



